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TECHNICAL  REPORT  SUMMARY 

Laser-induced  damage  to  initially  transparent  materials  represents 
a  serious  limitation  to  the  design  and  operation  of  high  power  lasers.  In 
order  to  assess  these  limitations  and  to  provide  standards  for  material 
performance,  it  is  necessary  to  obtain  reliable  measurements  of  damage 
intensities.  The  primary  goal  of  the  present  program  is  to  measure 
damage  intensities  in  various  transparent  solids  and  to  interpret  these 
results  ir.  terms  of  possible  damage  mechanisms.  Using  techniques  which 
we  developed  in  earlier  work,  we  have  been  able  to  isolate  and  to  study  the 
intrinsic  damage  process  of  electron  avalanche  breakdown.  Such  work 
establishes  measured  upper  limits  for  the  propagation  of  high  intensity 
light  in  solids  and  indicates  how  these  intrinsic  limits  change  with  laser 
pulse  and  material  characteristics.  We  have,  in  addition,  observed 
damage  from  small  absorbing  inclusions  and  have  begun  to  develop  tech¬ 
niques  which  can  be  applied  to  optical  materials  evaluation.  As  part  of 
the  present  effort,  a  C02  TEA  laser  was  designed  and  constructed  as  a 
TEMqq  mode,  10.6  p.m  source  for  optical  damage. 

A  study  of  intrinsic  damage  in  polycrystalline  and  disordered  ma¬ 
terials  has  been  completed  and  interpreted  in  terms  of  electron  avalanche 
breakdown.  By  comparing  the  optical  bulk  damage  fields  for  a  polycrystal, 
various  single-crystal  alloys,  and  two  glassy  solids  to  the  damage  fields 
for  the  respective  pure  single  crystals,  it  was  found  that  only  severe  lat¬ 
tice  disorder  such  as  present  in  completely  amorphous  fused  quartz  causes 
the  damage  fields  to  increase.  For  the  polycrystal  and  the  less  disordered 
systems,  the  damage  fields  are  the  same  as  those  of  the  pure  single  crystal. 
It  was  found,  in  addition,  that  high  concentrations  of  certain  types  of  atomic 
impurities  have  no  measurable  effect  on  intrinsic  damage  levels. 

Mease  -ements  have  been  made  of  the  pulse  width  independence  of 
intrinsic  damage  in  NaCl  at  1. 06  (im.  It  was  found  that  the  intrinsic 
damage  field  increased  by  nearly  an  order  of  magnitude  to  over  10 
volts/cm  as  the  laser  pulse  duration  was  decreased  from  10  ns  to  15  ps. 
These  results,  interpreted  in  terms  of  electron  avalanche  breakdown, 
are  the  first  measurements  of  intrinsic  damage  induced  by  mode-locked 

laser  pulses. 
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Optical  damage  to  dielectric  coatings  was  studied  with  weakly 
focused  ruby  laser  pulses  having  durations  between  20  ns  and  20  ps. 

Under  these  conditions  of  measurement,  damage  was  found  to  result 
from  highly  absorbing  inclusions.  Inclusions  with  diameters  ~  3.5  pm 
were  most  easily  damaged  with  20  ns  pulses  while  inclusions  with 
diameters  ~  0.4  pm  determined  damage  resistance  to  20  ps  pulses. 

Bulk  optical  damage  from  inclusions  was  observed  in  measure¬ 
ments  at  10.6  pm  on  chemical  vapor  deposition  grown  ZnSe.  The  in¬ 
clusion  damage  intensity  was  observed  to  be  statistical  in  nature  and, 
as  found  recently  in  studies  of  thin  film  damage,  varied  with  the  degree 
of  external  focusing.  A  statistical  model  of  inclusion  damage  was  used 
in  the  present  work  to  explain  these  observations  and  to  calculate,  from 
the  experimental  data,  an  average  inclusion  density.  The  intrinsic 
damage  intensity  for  ZnSe  was  also  determined. 

Finally,  a  classical  model  of  electron  avalanche  breakdown  was 
used  to  interpret  recent  experimental  data  on  intrinsic  laser  damage. 

The  model,  based  on  a  Fokker-Planck  equation  formula,  was  used  to 
show  the  relationship  between  dc  and  laser-induced  damage. 

The  present  work  on  intrinsic  damage  has  extended  over  knowledge 
of  electron  avalanche  breakdown.  It  may  now  be  possible  to  ref'ne  our 
models  of  this  damage  process  in  order  to  enable  us  to  accurately  predict 
the  effects  of  material  parameters  on  intrinsic  damage  fields  at  red  and 
near  infrared  frequencies.  Damage  measurements  should  be  extended 
to  higher  optical  frequencies  in  order  to  determine  if  damage  from  elec¬ 
tron  avalanche  breakdown  still  dominates  as  material  bandgaps  are  ap¬ 
proached  . 

Our  studies  of  inclusion  and  self-focusing  damage  can  serve  as  a 
model  for  future  materials  evaluation.  Since  inclusions  and  self-focusing 
normally  limit  the  intensity  of  light  that  materials  can  withstand,  such 
evaluation  is  essential  for  the  development  of  damage-resistant  optics 
and  for  the  design  of  damage-resistent  systems. 
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FOREWORD 


This  scientific  report  describes  work  performed  under  Contract 
No.  F  19623-73-C-0127  letween  1  January  1973  and  31  December  1973. 

The  report  was  assigned  a  Raytheon  internal  number  S-1660. 

Work  was  carried  out  at  the  Raytheon  Research  Division  in  Waltham, 
Massachusetts.  The  Principal  Investigator  until  June  30,  1973  was 
Dr.  Michael  Bass.  After  that  date  Dr.  David  Fradin  became  Principal 
Investigator.  Additional  experiments  were  performed  at  the  Naval  Research 
Laboratories,  Washington,  D.C.,  in  collaboration  with  Dr.  J.  P.  Letellier, 
and  at  the  Air  Force  Cambridge  Research  Laboratories,  Bedford,  Mass., 
in  collaboration  with  Dr.  David  Milam  and  R.  A.  Bradbury.  This  technical 
report  was  prepared  under  the  direction  of  Dr.  Fradin  in  collaboration  with 
D.  P.  Bua.  The  authors  have  benefitted  from  the  many  discussions  with 
Keimpe  Andringa,  Dr.  Frank  A.  Horrigan,  Dr.  Thomas  Deutsch  and 
Dr.  Robert  Rudko.  Miss  Carol  Christian  assisted  in  the  studies  of  Sec.  VII 
and  VIII.  Dr.  A.  Linz  generously  provided  several  of  the  sample  used  in  the 
damage  study  of  disorder  materials  and  J.  VanderSande  conducted  electron 
reflection  diffraction  measurements  for  that  study. 

This  report  was  submitted  by  the  authors  on  31  January  1974. 
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GENERAL  INTRODUCTION 


The  basic  objectives  of  this  program  was  to  analyze  data  obtained 
from  carefully  controlled  damage  measurements  at  laser  wavelengths 
between  0.69  and  10.6  urn  and  at  various  laser  pulse  durations.  Laser 
damage  thresholds  of  various  materials  of  interest  were  determined 
using  techniques  developed  during  the  past  year,  and  the  results  studied 
to  ascertain  the  mechanisms  responsible  for  damage.  Theoretical  studies 
of  intrinsic  laser  damage  from  avalanche  breakdown  were  conducted,  and 
a  new  technique  was  developed  to  measure  self-focusing  parameters  by 
inducing  optical  damage. 

The  effort  was  divided  into  six  major  areas:  They  were  (1)  design 
and  assembly  of  a  pulsed  CO^  TEA  laser  for  use  in  10.6  pm  damage  meas¬ 
urements,  (2)  a  study  of  the  effects  of  lattice  disorder,  polycrystallinity, 
alloying,  and  atomic  impurity  doping  on  intrinsic  damage  fields,  (3)  meas¬ 
urements  of  the  dependence  of  intrinsic  breakdown  fields  on  laser  pulse 
duration,  (4)  bulk  and  thin  film  optical  damage  from  inclusion  absorption, 

(5)  the  relationship  between  self -focusing  and  optical  damage,  and  (6) 
analytical  studies  of  laser-induced  electron  avalanche  breakdown. 

A  major  effort  in  the  present  program  was  the  design  and  construction 
of  a  COg  pulsed  laser.  By  using  Rogowski  shaped  electrodes  and  atmospheric 
resonator  pressures,  it  was  possible  to  produce  a  uniform,  large  volume 
discharge  inside  the  resonator  and  to  obtain  approximately  0.8  joules  of 
multimode  laser  output  with  pulse  durations  of  about  140  ns.  With  intra¬ 
cavity  aperturing  the  laser  intensity  distribution  was  restricted  to  a  gaussian 
mode  having  peak  powers  of  about  1  MW.  This  power  was  more  than  suffi¬ 
cient  for  optical  damage  studies. 

Intrinsic  optical  damage  was  investigated  in  a  cooperative  effort  in¬ 
volving  both  the  Raytheon  Research  Division  and  the  Naval  Research  Labora¬ 
tories.  The  work  at  Raytheon  investigated  the  possible  effects  of  selected  material 
parameters  on  intrinsic  damage  fields.  It  was  found  that  intrinsic  damage 
is  insensitive  to  atomic  impurities  and  lattice  disorder.  Only  the  severe 
lattice  disorder  that  is  found  in  a  truly  amorphous  solid  such  as  fused 
quartz  can  increase  intrinsic  breakdown  limits.  Alloyed  material  have 
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intrinsic  breakdown  intensities  intermediate  between  those  of  the  pure 
constituent.  These  results  are  shown  to  be  consistent  with  an  electron 
avalanche  damage  mechanism. 


Experiments  at  NRL,  in  collaboration  with  personnel  from  NRL 
and  Harvard  Lniversity,  probed  the  dependence  of  intrinsic  breakdown 
fields  on  laser  pulse  duration.  The  laser  intensity  necessary  to  induce 
damage  in  NaCT  was  found  to  increase  for  pulse  durations  •$.  1  ns.  This 
dependence  is  shown  to  be  consistent  with  dc  dielectric  breakdown  meas¬ 
urements  on  thin  samples  of  NaCl.  The  agreement  between  laser  and  dc 
results  underscores  the  fundamental  similarity  between  dc  and  laser- 
induced  electrical  breakdown  in  insulators. 

Studies  of  thin  film  damage  from  inclusion  absorption  were  con¬ 
ducted  by  the  Air  Force  Cambridge  Research  Laboratories  at  AFCRL’s 
Optical  Physics  Laboratory,  with  Raytheon  personnel  participating  in  the 
analysis  and  interpretation  of  results.  By  using  weakly  focused  ruby  laser 
pulses  of  varying  time  duration  and  by  carefully  studying  the  residual 
damage  morphology,  it  was  found  that  threshold  damage  resulted  from  in¬ 
clusions  whose  diameters  were  related  in  a  predictable  manner  to  laser 
pulse  width. 

Bulk  damage  from  inclusion  absorption  was  studied  by  using  the 
C02  TEA  laser  to  induce  damage  inside  chemical  vapor  deposition  grown 
ZnSe  polycrystals.  Consistent  with  previously  published  measurements 
of  thin  film  damage,  it  was  found  that  the  damage  intensity  varied  with 
the  degree  of  external  focusing.  A  simple  model  of  inclusion  damage 
was  developed  to  explain  this  dependence  and  is  used  to  infer  an  average 
inclusion  density  from  the  damage  data.  It  is  shown  that  inclusion  damage 
can  appear  to  be  statistical  in  character  because  of  the  random  placement 
of  inclusions  inside  the  sample. 

The  known  relationship  between  self-focusing  and  optical  damage 
at  low  laser  input  powers  is  used  as  the  basis  for  a  new  technique  to 
measure  self-focusing  parameters.  As  an  illustration  of  the  use  of  the 
technique,  the  self-focusing  parameters  for  sapphire  are  determined. 
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finally,  classical  model  of  avalanche  breakdown  is  presented  in 
order  to  provide  a  semiquantitative  framework  for  interpreting  the  in¬ 
trinsic  damage  studies.  This  model,  based  on  a  Fokker-Planck  equation 
formalism,  is  an  extension  of  earlier  analysis  at  Raytheon. 


II. 


DESIGN  AND  CONSTRUCTION  OF  A  CO,  TEA  LASER 


A .  Introduction 

A  major  objective  of  the  present  program  was  the  design  and 
construction  of  a  CC>2  transverse-excited-atmospheric  (TEA)  laser. 

This  work  was  divided  into  three  phases: 

1.  selection  of  a  gas  discharge  design 

2.  construction  of  an  operating  laser  resonator 

3.  application  of  design  refinements  of  the  resonator  optics 

to  produce  TEMqq  output  with  sufficient  power  and  spectral 
purity  to  perform  optical  damage  measurements. 

The  three  phases  have  been  completed  and  the  laser  has  been  de¬ 
livered  to  the  Air  Force  Cambridge  Research  Laboratories.  A  double 
Rogowski  electrode  configuration1  was  selected  for  our  use.  Using  a 
published  electrical  design  with  minor  modifications,  we  were  able  to 
achieve  a  stable,  uniform  discharge,  and  to  produce  approximately 
0.2  joules  of  TEMqo  mode  output  at  10.6  |xm.  Single-longitudinal  mode 
operation  was  also  achieved  by  using  an  intracavity  germanium  etalon. 

B.  Resonator  Design 

High  peak  powers  can  be  obtained  with  CC>2  lasers  by  use  of  atmos¬ 
pheric  pressure  and  fast  electric  discharges.3  There  are  two  reasons  for 
using  transverse  excitation  at  high  pressures  instead  of  a  longitudinal  dis¬ 
charge.  First,  lower  voltages  are  required  to  produce  gas  breakdown 
since  the  discharge  path  is  reduced.  In  addition,  the  discharge  impedance 
is  lower  for  transverse  excitation  so  that  energy  can  be  rapidly  injected 
into  the  discharge  volume.  This  second  advantage  is  important  because 
it  allows  the  laser  to  be  excited  in  times  short  compared  to  the  excited - 
state  lifetime  of  the  C02  molecule,  a  lifetime  equal  to  about  10  fis  at 
atmospheric  pressure.4  As  a  result  of  this  rapid  excitation,  the  laser 
can  be  Q-switched  without  using  such  intracavity  devices  as  electro-optic 
shutters  or  rotating  mirrors. 
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A  number  of  schemes  for  transverse  excitation  have  been  discussed 
in  the  literature.  *  ’  Two  of  them,  pin  electrodes  and  Rogowski  elec¬ 
trodes  with  trigger  wires,  were  compared  in  our  laboratory.  The  latter 
design  was  chosen  because  il  offered  the  advantages  of  a  large  discharge 
volume,  comparative  simplicity,  and  a  relatively  low  level  of  radiated  rf 
noise. 


The  laser  was  constructed  with  two  parallel,  graphite  electrodes 

n 

shaped  to  approximate  a  Rogowski  profile.  Two  thin  (9-mil)  tungsten 
trigger  wires  running  parallel  to  the  electrodes  and  placed  just  outside 
the  lasing  region  (Fig.  1)  provided  preionization  of  the  discharge  volume* 
with  the  main  discharge  occurring  across  the  graphite  electrodes. 
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Chang  and  Wood  have  published  a  simple  electrical  arrangement 
for  this  type  of  laser,  and  we  adopted  their  design  with  minor  modifica¬ 
tions.  Figure  2  shows  a  schematic  of  this  design  with  the  circuit  clement 
values  that  we  used  in  our  system.  To  produce  the  discharge,  capacitors 
Cj  and  C2  were  initially  charged  to  about  30  kV.  A  short  voltage  spike 
from  a  trigger  transformer  fired  spark  gap  SG,  causing  current  ij  to 
flow.  The  spark  gap  was  a  commercial  device  (Tachisto)  especially 
modified  for  our  use.  Because  of  the  spark  gap  inductance,  ij  started 
to  oscillate  and  drive  the  trigger  wires  T  and  cathode  K  to  negative 
voltages .  Wi  hen  the  voltage  at  T  began  to  drop  below  ground,  an  arc 
discharge  was  produced  between  T  and  the  resonator  anode  A.  UV  light 
from  this  arc  caused  preionization  in  the  gas  volume  between  the  laser 
electrodes  K  and  A  and  helped  to  initiate  a  large-volume  uniform  main 
discharge.  This  discharge  rapidly  damped  out  the  ringing  current  and 
absorbed  most  of  the  energy  stored  in  the  circuit.2 

The  discharge  characteristics  were  influenced  by  the  gases  flowing 
in  the  resonator.  We  used  approximately  800  torr  of  N9:CO„:He  mixed 

“  u 

in  the  ratio  3.2:3.2:10.5.  This  mixture  produced  both  a  uniform  glow 

disrhrage  and  optimized  laser  output.  No  vacuum  pump  was  used  in  this 
2 

system.  Instead,  the  gases  are  passed  through  an  exhaust  in  the  reson¬ 
ator  box  at  a  total  rate  of  about  33  liters/minute. 
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Anode 


Cathode 


Trigger  wires 


Rogowski  Shaped  Graphite  Electrodes.  Each  electrode  was 
polished  to  remove  sharp  edges.  The  trigger  wires  were 
9  mil  tungsten. 
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The  physical  design  of  the  resonator  was  also  modeled  after  Ref.  2. 
Figure  3  is  a  sketch  of  the  resonator  showing  the  placement  of  the  charging 
capacitors  and  spark  gap  and  the  design  of  the  electrode  mounts.  A  photo¬ 
graph  of  the  actual  laser  is  reproduced  in  Figs.  4  and  5.  Because  of  the 
use  of  thick  resonator  walls,  the  laser  could  be  operated  with  other  gases 
at  well  below  atmospheric  pressure. 

Initial  tests  on  this  resonator  (Fig.  6)  produced  about  0.7  joules  of 
multimode  laser  output  at  10.6  urn  with  peak  powers  of  about  0.6  mW  for 
an  input  electrical  energy  of  29  joules.  By  optimizing  the  mirrors  (ger¬ 
manium  mirror  100  percent  reflective  at  10.6  urn,  10  meter  radius,  and 
a  germanium  mirror  80  percent  reflective  at  10.6  nm,  flat)  we  were  able 
to  obtain  0.8  joules  of  multimode  energy  with  a  cavity  spacing  of  1.  2  meters 
and  an  input  electrical  energy  of  26.5  joules.  Firing  at  about  1  pps,  the 
pulse  energy  stability  was  about  5  percent.  (Figure  7) 

C .  Optical  Design  of  the  Resonator 


Performance  requirements  for  the  CO2  laser  were  determined  by 
its  intended  use  as  a  light  source  for  optical  damage  experiments.  It  was 
essential  for  our  purposes  that  the  laser  operate  in  a  TEMQQ  mode  and 
that  it  achieve  peak  powers  in  excess  of  about  150  kW.  It  was  desirable 

that  the  laser  output  contain  a  single  longitudinal  mode  or  that  it  be  made 
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completely  time-resolved  as  in  the  work  of  Yablonovitch,  particularly  if  the 
damage  mechanism  involves  inclusion  absorption^  with  nanosecond  or  sub¬ 
nanosecond  thermal  time  constants.  Spectral  purity  is  less  critical  for  elec¬ 
tron  avalanche  breakdown^  as  the  effect  of  high-frequency  mode  beating  tends 
to  average  out. 

A  number  of  optical  configurations  can  be  used  to  produce  TEMqq 
mode  operation.  The  simplest  is  a  stable  resonator  having  mirrors  with 
radii  of  curvature  larger  than  the  resonator  length  and  having  an  intra¬ 
cavity  aperture  constraining  the  cavity  Fresnel  number. 

Another  resonator  configuration  we  considered  was  an  unstable 
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resonator  with  annular  output  coupling.  Initial  tests,  however,  clearly 
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indicated  that  the  first  method,  using  a  curved  mirror  and  an  aperture 
in  a  stable  resonator,  would  supply  us  with  more  than  enough  energy 
for  our  experiments.  The  beam  profile  obtained  with  a  stable  resonator 
was  more  conducive  to  damage  work  since  the  output  intensity  distribu¬ 
tion  of  the  TEMqq  stable  resonator  was  verv  easily  handled  analytically. 
Since  our  earlier  damage  work  had  been  conducted  with  Gaussian  beams, 
comparative  analytical  interpretations  were  clearer  with  the  stable  re¬ 
sonator  geometry. 

A  third  resonator  configuration  considered  was  a  stable  resonator 
with  intracavity  optics  to  enlarge  the  mode.  This  technique  was  compli¬ 
cated  and  unnecessary  for  our  purposes. 

Figure  B  show’s  the  laser  cavity  with  aperture  installed.  As  the 
aperture  size  was  decreased,  the  pulse  time-structure  (as  observed  with 
a  GHz  pyroelectric  detector)  was  seen  to  improve  to  the  point  where  the 
laser  output  became  simultaneously  mode-locked  (Fig.  9).  This  change 
in  time-structure  was  a  good  indication  that  the  laser  was  operating  TEMqq 
mode.  The  fact  that  the  beam  intensity  distribution  is  gaussian  was  demon¬ 
strated  by  sweeping  a  pinhole  across  the  beam  diameter  while  monitoring 
the  energy  transmitted  through  the  pinhole.  A  plot  of  the  log  of  transmitted 
energy  versus  the  beam  radius  squared  produced  a  straight  line  (Figs.  10 
and  11).  With  an  input  of  electrical  energy  of  26  joules,  an  output  of 
120  mJ  of  TEMqq  mode  laser  light  was  obtained. 

Even  though  longitudinal  mode  selection  was  not  necessary  for  our 
damage  studies,  some  short  experiments  were  conducted  with  this  laser 
to  determine  if  longitudinal  mcde  selection  could  be  easily  accomplished. 

With  the  laser  operating  TEMqq  mode,  longitudinal  mode  selection 

was  accomplished  by  using  an  intracavity  germanium  etalon  aligned  normal 
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to  the  cavity  axis  '  (Fig.  12).  This  same  technique  had  been  used 

previously  to  restrict  the  number  of  lasing  modes  in  ruby  and  Nd:YAG 

lasers.  In  a  solid  state  laser  an  intracavity  etalon  can  be  formed  by  the 

highly  reflecting  end  mirror  and  the  plane,  antireflection-coated  front 
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surface  of  the  laser  rod.  The  high  gain  of  the  medium  enhances  the 
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effective  reflectivity  of  the  rod  face  sufficiently  to  produce  an  effective 

etalon  with  a  finesse  of  about  6  and  a  free  spectral  range  comparable  to 

that  of  the  total  resonator.  In  this  manner,  the  number  of  longitudinal 

modes  in  YAG  and  ruby  lasers  were  restricted  to  four  or  less.  In  COg 

TEA  lasers,  an  uncoated  germanium  flat  positioned  between  the  discharge 
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region  and  the  output  mirror  produces  similar  spectral  selection. 

Figure  13  is  a  picture  of  the  laser  output  obtained  with  the  germanium 
etalon  in  place.  The  photograph  was  made  with  a  500  MHz  pyroelectric 
detector /oscilloscope  combination. 

While  the  germanium  etalon  did  produce  single  frequency  output, 
it  also  caused  a  loss  in  the  laser  output  energy  and  output  stability.  The 
loss  in  output  stability  was  apparently  mechanical  in  origin  and  could  have 
been  solved  by  mounting  the  etalon  in  a  more  stable  mount.  In  order  to 
improve  the  resonator  stability,  different  materials  with  lower  indices 
of  refraction  and  hence  lower  surface  reflectivity  were  tested.  NaCl  did, 
in  fact,  lead  to  greater  laser  stability,  but  the  reflectivity  from  its  sur¬ 
faces  was  so  low  that  single  frequency  operation  could  not  be  attained. 

ZnSe,  whose  index  of  refraction  lies  between  that  of  sodium  chlorid 
and  germanium,  appeared  to  present  a  good  compromise  between  spectral 
selectivity  and  laser  stability.  We  did  not  adequately  pursue  the  use  of 
ZnSe,  however,  because  the  optical  quality  of  the  material  available  was 
low.  Significant  pulse  smoothing  was  nonetheless  observed  with  the  ZnSe 
etalon.  Experiments  with  SFp  and  chlorotrifluororethylene  gases  were 
not  successful  because  these  gases  produced  only  selected  line  absorption 
without  an  intracavity  grating:  their  use  merely  forced  the  laser  to  change 
its  oscillation  to  a  different  vibrational-rotational  line. 

Pulsed  CC>2  lasers  may  operate  on  more  than  one  rotational  line. 
Since  the  lines  of  the  P  branches  of  either  rotational  band  are  separated 
by  about  50  GHz,  simultaneous  lasing  on  more  than  one  rotational  line 
will  produce  ultrahigh  frequency  beating.  Such  fast  time  structure  cannot 
affect  the  results  of  optical  damage  measurements.  For  this  reason,  no 
special  attempt  was  made  to  restrict  lasing  to  a  single  rotational  branch. 


D.  Conclusion 


"  8  haVe  COmpleted  the  design  and  construction  of  an  operating  CO 
TEA  laser  that  more  than  meets  the  requirements  needed  to  perform  op-  2 
deal  damage  experiments.  Figure  14  shows  the  arrangement  we  used  in 
our  expertments.  The  laser  operated  TEM  mode  with  approximately 
120  mJ  of  energy  of  which  60  mj  was  contained  in  the  main  pulse  and  ap¬ 
proximately  60  mJ  in  the  tail.  A  combination  of  a  rotatable  stack  polarizer 
(two  germanium  flats)  for  fine  adjustment  and  two-inch  square  pieces  of 
calcium  fluoride  attenuators  for  coarse  adjustments  were  used  to  control 
the  laser  power  incident  on  test  sample.  The  output  beam  was  focused 
w.th  antireflection  coated  germanium  lenses  having  various  focal  lengths. 
All  energy  measurements  were  made  with  a  fast  response  Jen/Tec  energy 
meter.  Pulse  response  was  monitored  with  a  500  pic  second  pyroelectric 
detector  manufactured  by  Molectron  (model  P5E)  and  a  Tektronix  type 
7904  oscilloscope. 
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111  •  effects  of  lattice  disorder  on  the  intrinsic 

OPTICAL  DAMAGE  FIELDS  OF  SOLIDS 


A .  Introduction 

Intrinsic  damage  fields  have  previously  been  determined  only  for 
pure  single  crystals.  Recent  interest  has  developed,  however,  in  dis¬ 
ordered  solids  as  possible  new  materials  for  high  power  laser  optics. 
Disordered  solids,  particularly  polycrystals  and  alloys,  are  attractive 
because  of  their  superior  physical  properties.  In  addition,  certain  ma¬ 
terials  such  as  ZnSe  are  more  easily  produced  in  polycrystalline  form 
than  as  single  crystals.  There  are,  then,  practical  reasons  for  meas¬ 
uring  the  intrinsic  breakdown  fields  of  disordered  solids. 

A  more  fundamental  reason  for  such  measurements  concern  the 
nature  of  the  intrinsic  damage  process.  It  has  been  demonstrated  that 
electron  avalanche  breakdown  is  responsible  for  intrinsic  damage  of 
transparent  solids.  Since  our  previous  work  had  established  that  ef¬ 
fective  electron  collision  times  are  on  the  order  of  10" 16  to  lO-15  sec, 
it  appears  that  the  hot  electrons  involved  in  the  avalanche  are  changing 
their  moments  rapidly  by  effectively  colliding  with  every  atom  in  their 
paths.  This  fact  implies  that  only  severe  lattice  disorder  should  affect 
the  avalanche.  It  is  argued  in  Sec.  B,  in  fact,  that  the  intrinsic  damage 
field  does  not  change  unless  disorder  exists  on  the  scale  of  10  atomic 
units  or  less.  If  such  severe  disorder  exists,  it  will  be  more  difficult 
to  heat  the  electron  population  and  hence  the  damage  field  will  rise.  A 
measurement  of  the  intrinsic  damage  in  disordered  materials,  therefore, 
is  a  test  of  our  physical  conception  of  electron  avalanche  breakdown. 

We  have  measured  the  intrinsic  optical  damage  fields  of  three 
systems  having  different  degrees  of  disorder  -  polycrystal  KC1,  single¬ 
crystal  KCl-KBr  alloys,  and  glass  solids.  It  was  found  that  extreme 
lattice  disorder  such  as  present  in  a  highly  amorphous  material  such 
as  fused  quartz  causes  the  damage  field  to  increase  whereas  less  severe 
disorder  does  not  measurably  affect  the  breakdown  strength.  This  result 
is  consistent  with  our  predictions  based  on  avalanche  breakdown  as  being 
the  damage  mechanism.  It  demonstrates,  in  addition,  that  the  intrinsic 
limits  for  light  propagation  in  single  crystals  and  polycrystals  of  the 
same  materials  are  identical. 
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B*  Effects  of  Lattice  Disorder  on  the  Intrinsic  Optical 
Damage  Fields  of  Solids 

Measurements  of  the  effects  of  lattice  disorder  on  the  intrinsic 
laser-induced  damage  fields  of  transparent  solids  are  reported.  Ex¬ 
treme  lattice  disorder  such  as  present  in  a  highly  amorphous  material 
such  as  fused  quartz  causes  the  damage  field  to  increase  whereas  less 
severe  disorder  does  not  measurably  affect  the  breakdown  strength. 

This  is  consistent  with  an  electron  avalanche  intrinsic  damage  mechan¬ 
ism. 


Studies  of  dc  dielectric  breakdown  have  shown  that  severe  lattice 
disorder  can  increase  the  electric  strength  of  solids. 16  Since  it  has 
recently  been  observed  that  intrinsic  optical  damage  in  transparent 
solids  appears  to  develop  from  the  same  mechanism  that  is  responsible 
for  the  dc  dielectric  breakdown,  9*  10* 17  severe  lattice  disorder  should 
also  increase  the  optical  damage  fields  of  solids.9 

Measurements  are  reported  in  this  paper  of  optical  bulk  damage 
in  three  systems  having  different  degrees  of  disorder  --  polycrystal  KC1, 
single-crystal  KBr-KCl  alloys,  and  glassy  solids.  In  each  case  the  dam¬ 
age  field  for  the  disordered  system  was  compared  to  the  optical  strength 
of  the  corresponding  crystal.  It  was  found  that  severe  lattice  disorder 
such  as  present  in  the  completely  amorphous  system  (fused  quartz)  causes 
the  intrinsic  optical  damage  fields  to  increase.  This  observation  is  con¬ 
sistent  with  dc  breakdown  experiments 16  and  with  simple  models  of  elec¬ 
tron  avalanche  breakdown.  19 

The  laser  system  and  the  experimental  techniques  used  here  are 
described  in  Ref.  10.  A  Q-switched,  TEMQQ  Nd:YAG  laser  was  used  to 
induce  damage  inside  the  bulk  of  various  samples.  Self-focusing  was 
avoided  by  restricting  the  laser  input  power,  9» 10  and  damage  from  in¬ 
clusions  was  distinguished  from  intrinsic  damage  by  examining  both  the 
morphology  of  the  damage  sites  '  9  and  the  temporal  shape  of  light 
pul,  es  transmitted  through  the  sample.  Only  data  obtained  from  in¬ 
trinsic  damage  events  were  considered  in  the  present  work.  The  damage 
field  was  defined  as  that  value  of  on-axis,  root-mean-square  electric  field 
necessary  to  induce  damage  on  a  single  shot  with  a  probability  of  0.5.10»20 
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The  damage  field  of  the  large-grain  (20  nm)  polycrystal  uas  the 
same  as  that  measured  in  the  single  crystal,  and  the  damage  fields  for 
the  alloys  were  intermediate  between  the  damage  fields  of  the  constituents. 
In  quartz,  on  the  other  hand,  the  disordered  form  was  noticeably 
stronger  than  the  crystal,  the  ratio  of  damage  intensities  being  5  ±  1. 

This  ratio  is  identical  to  the  corresponding  ratio  of  surface  damage  fields 
measured  in  previous  work. 


Damage  data  were  taken  on  a  second  amorphous -crystal  pair, 
Ba2^S^e2®7  (BMGO).  Several  measurements  had  indicated  that  the 
BMGO  glass  was  distorted.  An  X-ray  Laue  pattern  of  the  glass  gave 
no  evidence  of  crystallinity.  Similarly,  the  directions  of  the  stress- 
induced  fractures  that  normally  accompany  optical  damage^  (Fig.  15) 
were  random  in  the  amorphous  BMGO  as  expected  for  a  glass,  where¬ 
as  fracture  in  the  crystal  occurred  along  well-defined  crystal  directions. 

Disorder  on  a  local  scale  in  BMGO  glass  was  evident  in  the  lumines- 

21 

cence  studies  of  Munasinghe  and  Linz.  When  our  sample  was  carefully 
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studied  with  electron  reflection  diffraction,  however,  it  was  found  that 
regions  ( ^.  2000A)  of  crystal  ordering  existed.  Within  the  volume  ir¬ 
radiated  by  the  laser,  therefore,  the  BMGO  glass  was  not  completely 
amorphous.  It  was  actually  a  mixed  crystal-amorphous  phase.  For 
the  BMGO  system,  the  damage  fields  of  the  glass  and  the  crystal  were 
the  same  within  experimental  error. 


Table  I  and  Fig.  16  summarize  the  measured  optical  damage  fields. 

It  is  to  be  expected  that  the  large-grain  polycrystal  should  have 
the  same  damage  field  as  the  single  crystal.  The  average  grain 
(20  pm)  and  the  laser  focal  volume  were  comparable  so  that  in  the  high 
intensity  region  near  the  beam  axis  where  breakdown  is  observed  to  initiate, 
the  sample  looks  like  a  single  crystal. 

By  a  simplified  argument  we  can  estimate  the  degree  of  disorder 

necessary  to  affect  the  breakdown  strength.  Classical  theories  of  aval- 
1 9 

anche  breakdown  predict  that  the  dynamics  of  electrons  with  energies 
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greater  than  the  longitudinal  optical  (LO)  energy  determine  the  charac¬ 
teristics  of  the  avalanche.  The  LO  energy  in  the  alkali  halides  corres¬ 
ponds  to  electron  momenta,  k,  of  about  0.  1  times  the  reciprocal  lattice 
vector,  G.  Thus  the  important  electrons  have  k  ^  0. 1G  and  their  inter¬ 
action  with  the  lattice  will  be  dominated  by  phonons  having  wave  number 
>  0.  1G.  Because  such  lattice  vibrations  have  wavelengths  equal  to  10 
lattice  constants  or  less,  this  simple  model  suggests  that  the  damage 
field  should  not  be  affected  by  disorder  unless  the  disorder  appears  on 
the  scale  of  about  10  lattice  constants  (~  50A).a 

Amorphous  systems  may  be  disordered  on  such  a  scale.  Our 
observation  that  fused  quartz  is  more  resistant  to  damage  than  crystalline 
quartz  is,  therefore,  consistent  with  the  argument  just  presented.  The 
data  on  the  Ba2MgGe20^  system  is  also  explained  by  this  argument  since 
that  system  had  regions  within  the  focal  volume  which  were  locally  ordered 
on  the  scale  of  >  >  10  lattice  constants.  Electron  avalanche  is  apparently 
initiated  in  these  crystalline  regions. 

The  breakdown  fields  for  various  single  crystal  alloys  of  KBr 
and  KC1  are  summarized  in  Fig.  16.  It  is  seen  that  the  damage  strengths 
of  the  alloys  are  intermediate  between  those  of  the  constituents.  There  is 
no  evidence  that  alloying  caused  disorder  S'lfficient  to  increase  the  break¬ 
down  strength  of  these  materials.  The  variation  of  damage  strength  with 
composition  can  be  qualitatively  understood  by  noting  that  many  material 
parameters  such  as  bandgap,  lattice  constant,  dielectric  constant,  and 
phonon  frequencies  have  values  intermediate  between  those  of  the  con¬ 
stituents?'^  According  to  simple  models  of  avalanche  breakdown,1®  the 
breakdown  field  depends  on  these  various  material  parameters  so  that  it 
is  reasonable  to  expect  the  breakdown  strengths  of  the  alloys  to  also  be 
intermediate  between  those  of  the  constituents. 


aThe  larger  breakdown  field  for  the  amorphous  quartz  can  be  understood 
by  another  physically  equivalent  argument.  As  noted  in  Refs.  9  and  10, 
the  rate  of  energy  input  into  the  electron  population  decreases  with  de¬ 
creasing  electron  mobility  u.  In  low  fields  n  is  smaller  in  highly  dis¬ 
ordered  systems  than  it  is  in  crystals  with  the  same  chemical  composition. 
If  the  effective  mobility  in  disordered  media  is  also  lower  in  the  ultra- 
high  fields  characteristic  of  avalanche  breakdown,  then  it  should  be  more 
difficult  to  heat  the  electron  distribution  in  fused  quartz,  and,  as  ob¬ 
served,  fused  quartz  should  be  more  damage  resistant  than  crystalline 
quartz. 
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It  thus  appears  that  extreme  lattice  disorder  such  as  present  in 
highly  amorphous  systems  has  a  measurable  effect  on  intrinsic  damage 
fields.  This  results  is  qualitatively  predicted  by  simple  models  of 
electron  avalanche  breakdown. 

The  advice  and  assistance  of  D.  Bua,  R.  Newberg,  D.W.  Howe, 

C.  Willingham  and  R.R.  Monchamp  are  gratefully  acknowledged.  A.  Linz 
generously  provided  several  of  our  samples  and  J.  VanderSande  conducted 
the  electron  reflection  diffraction  studies. 
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IV. 


DEPENDENCE  OF  LASER  INDUCED  BREAKDOWN  FIELD 
STRENGTH  ON  PULSE  DURATION 


A .  Introduction 

Intrinsic  damage  fields  from  avalanche  breakdowns  have  been  deter¬ 
mined  for  a  number  of  materials  using  Q-switched  laser  pulses.  9’  10,  24  It 
was  found  that  the  damage  fields  are  nearly  independent  of  pulse  duration 
tp  when  tp  >  5  ns.  When  high  frequency  (>  GHz)  mode  beating  occurs, 
however,  the  resulting  time  structure  is  averaged  out  by  the  damage 
process,  thus  suggesting  that  characteristic  development  times  of  the 
electron  avalanche  are  on  the  order  of  10-9  sec. 

Measurements  of  avalanche  breakdown  with  subnanosecond  laser 
pulses  have  not  been  made  prior  to  the  present  work.  Since  laser  break¬ 
down  and  dc  dielectric  breakdown  appear  to  result  from  the  same  funda¬ 
mental  process,  data  at  dc  can,  in  principal,  be  used  to  estimate  the  pulse 
duration  dependence  to  the  laser  induced  avalanche.  Yablonovitch  and 
Bloembergen  have,  in  fact,  made  such  estimates  based  on  measured 
values  of  dc  breakdown  for  thin  samples  of  NaCl.  Their  estimates,  how¬ 
ever,  were  based  on  uncertain  limits  of  an  effective,  high-field  electron 
drift  velocity  and  ignored  important  experimental  uncertainties  such  as 
space  charge  development  which  can  influence  dc  breakdown  fields  but 
not  laser  breakdown  fields. 

There  is  no  reliable  data  in  the  literature,  therefore,  on  the  laser 
pulse  width  dependence  of  intrinsic  optical  damage.  Such  data  is  useful 
to  our  understanding  of  avalanche  breakdown  because  it  provides  infor¬ 
mation  on  the  time  development  of  the  electron  avalanche.  Such  data  also 
has  practical  importance  in  that  it  establishes  measured  upper  limits  for 
the  propagation  of  mode  locked  laser  pulses  in  solids. 

We  have  measured  the  intrinsic  breakdown  strength  of  NaCl  at 
1.06  pm  using  four  pulse  durations  between  10.3  ns  and  15  ps.  It  was 
found  that  the  intrinsic  breakdown  field  increases  by  a  factor  of  about 
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6  to  over  107  V/cm  as  the  pulse  duration  is  changed  by  about  three 
orders  of  magnitude.  From  these  measured  damage  fields,  an  ioni¬ 
zation  rate  for  the  electron  avalanche  can  be  inferred  and  compared 
to  the  predictions  of  Yablonovitch  and  Bloembergen. 


i 
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B.  DEPENDENCE  OF  LASER  INDUCED  BREAKDOWN  FIELD 
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and 
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Abstract 


Field  strengths  at  which  optical  damage  is  initiated  in  NaCl  have 
been  measured  with  a  mode-locked  Nd: YAG  laser  with  pulse  durations  of 
15  and  300  picoseconds.  Comparison  with  previously  reported  data  with 
a  Q-switched  laser  shows  that  the  field  strength  required  for  intrinsic 
optical  damage  increases  by  almost  one  order  of  magnitude  from  106  v/cm 
at  10  sec  to  over  10  v/cm  at  1.  5  X  10-11  sec.  This  is  in  qualitative 

agreement  with  published  estimates  based  on  the  electron  avalanche  break¬ 
down  mechanism. 
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Damage  produced  by  laser  beams  of  high  intensity  in  transparent 
materials  has  been  studied  intensively  for  many  years.  It  is  only  re¬ 
cently,  however,  that  the  effects  of  absorbing  inclusions  and  self-focusing 
have  been  carefully  eliminated,  and  the  intrinsic  breakdown  mechanism 
in  transparent  condensed  dielectric  media  has  been  positively  identified 
as  electron  avalanche  ionization.9"10,18'20  On  the  basis  of  this  mechan¬ 
ism  and  known  characteristics  of  dc  breakdown  by  avalanche  ionization, 
Yablonovitch  and  Bloembergen25  predicted  a  characteristic  dependence 
of  breakdown  field  strength  on  laser  pulse  duration.  In  this  note  an  experi¬ 
mental  determination  of  this  dependence  in  the  subnanosecond  regime  is 
presented,  which  turns  out  to  be  in  general  agreement  with  those  predictions. 


The  measurements  were  performed  by  focusing  mode-locked  YAG: 

Xd  laser  pulses  having  durations  of  15  and  300  picoseconds  inside  a  single 

crystal  of  NaCl.  Because  the  experimental  procedures  used  in  the  present 

work  were  identical  to  those  used  in  recent  studies  with  a  Q-switched  YAG: 

Xd  laser,  the  subnanosecond  measurements  can  be  directly  compared  to 

the  results  of  these  studies.  It  is  found  that  the  intrinsic  breakdown  field 

increased  by  almost  an  order  of  magnitude  to  over  107  volts/  cm  as  the  laser 

pulsewidth  was  decreased  from  10  nanoseconds  to  15  picoseconds.  Previous 

data  with  pulse  widths  typical  of  Q-switched  lasers  had  demonstrated  that 

the  dependence  on  pulse  width  in  the  nanosecond  regime  is  small  so  that  rela- 

tively  little  information  on  the  time  dependence  of  the  avalanche  process  was 

available.  Our  results  demonstrate  again  that  lasers  can  be  used  to  measure 

properties  of  dielectric  breakdown  which  are  difficult  to  obtain  by  dc  tech- 
26-28 

mques. 


A  pulsewidth  dependence  to  optical  damage  in  thin  films  was  first  observed 

Idoc  cS*  Bli®s1and  D*  Milam»  4th  ASTM  Symp.  on  Damage  in  Laser  Materials, 
MBS  Spec.  Pub.  372,  108  (1972). 

The  mechanism  of  damage  was  not  identified,  however.  Since  it  is  known 
that  thin  films  can  have  a  large  residual  absorption  and  significant  structural 
defects  and  that  both  effects  can  change  the  damage  characteristics  of  sur¬ 
faces,  it  is  not  clear  that  the  change  in  damage  fields  which  were  observed 
bv  Rliss  and  Milam  reflect  a  property  of  an  intrinsic  bulk  damage  mechanism. 
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field  was  well  defined  and  did  not  change  as  different  regions  of  the 
sample  were  probed  and  lenses  with  different  focal  lengths  were  used 
to  focus  the  radiation.  Also,  only  one  spark  occurred  with  each  dam¬ 
aging  laser  pulse,  and  the  spark  always  appeared  to  form  at  the  geom¬ 
etrical  focal  plane.  These  observations  contrasted  with  those  obtained 
under  conditions  whore  inclusion  damage  was  seen.30  it  was  therefore 
concluded  that  except  possibly  for  occasional  damage  sites,  inclusion 
damage  was  absent  in  highly  pure  NaCl  under  the  conditions  of  our  mea¬ 
surement  as  it  had  been  in  previous  work10where  Q-switched  lasers  were 
used. 


Tablo  II  summarizes  the  results  of  the  present  measurements  and 
those  of  Ref.  10.  An  increase  in  breakdown  strength  was  observed  as  the 
duration  of  the  laser  pulse  was  decreased.  As  the  pulse  duration  was 
changed  from  10.3  ns  to  15  ps,  there  was  a  total  change  by  a  factor  of 
5.8  in  damage  field  strengths  or  a  factor  of  33  in  damage  intensity.  The 

experimental  points  are  plotted  in  Fig.  17,  together  with  some  semi-empirical 
predicted  curves  taken  from  Ref.  2  5. 

In  the  subnanosecond  regime  electron  diffusion  and  self-trapping 

m  .y  he  ignored.  The  density  of  the  conduction  electrons  in  the  avalanche 
changes  exponentially  with  time, 


N(t)  =  NQ  exp[J  a  (E)dt]  =  NQM  (t). 


(1) 


When  the  electron  density  exceeds  about  10 18  cm'3,  requiring  a  multipli¬ 
cation  factor  of  roughly  Mc  ~108,  breakdown  is  said  to  occur.  According 
to  Eq.  (1)  the  ionization  rate  o-(E) 


follows: 


rms 


is  related  to  the  pulse  duration  as 


o  (E  ) 
rms 


■1 


t  {  n  M 


18  t 


-1 


(2) 


This  relation  has  been  used  to  convert  the  quantity  o(E)  used  along  the 
vertical  axis  in  the  figure  of  Ref.  25  to  our  figure  which  uses  t  '  1 .  We  have 
shifted  the  curves  along  the  horizontal  axis  to  obtain  agreement  with  the 
experimental  values  for  the  breakdown  field  Erms  for  the  long  pulses.  It 
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RMS  OPTICAL  ELECTRIC  FIELD  (MV/cm) 


Sh°l,ld  b°  n°,ed  ,ha'  lh0  dr  '>rc,kdow„  Holds  quoted  in  Refs.  2G-28,  on 
vhioh  the  curves  in  Ref.  25  were  based,  are  about  a  factor  2.3  lower 

la"  °UI  laSLr  lalui"’  f»f  Erms-  There  are  reasons  to  suspect  that  this 
is  due  to  systematic  errors  in  the  dc  experiment.  First,  this 
lac  lor  or  2.3  is  nearly  the  same  for  all  nine  alkali  halides  studied  at 
10.  f,  um  and  1.  Of,  um,  and  second,  dc  field  values  are  average  values 
Without  regard  for  field  inhomogeneities  from  space  charge  effects  which 
can  be  important  in  dc  experimental *  31  - 32  Consistent  with  the  approach 
of  Ref.  9  and  10,  it  is  therefore  more  meaningful  to  compare  trends  in 
breakdown  fields  as  different  materials  are  investigated  or  as  parameters 
are  varied  than  to  compare  the  absolute  values  of  breakdown  fields. 

It  is  seen  that  the  experimental  points  fall  close  to  the  upper  curve 

°'  l!°f'  25,  Which  was  d<fived  on  the  assumption  that  the  mobility  in  the 
hoi  electron  gas  is  independent  of  Ermt..  Quantitative  agreement  should 
no.  be  emphasised ,  however,  because  the  analysis  is  based  on  assumptions 
wh,ch  are  questionable  over  the  range  of  damage  fields  considered.  As 
mentioned  before,  space  charge  and  electrode  effect  may  considerably  al¬ 
ter  the  interpretation  of  the  dc  results19  and  consequently  to  transition 
from  b.dc  to  Erms  along  the  horizontal  axis.  Furthermore  the  breakdown 
field  at  the  shortest  pulse  durations  becomes  so  high  that  it  is  possible 
that  frequency-dependent  tunneling  (or  multiphoton  ionization)  takes  over 
as  an  intrinsic  damage  mechanism.  In  Ref.  25  it  was  estimated  that  this 
change  over  would  occur  for  fields  higher  than  2  X  If,7  v/cm  or  pulse 
durations  shorter  than  one  picosecond.  If  this  estimate  is  inaccurate  and 
this  other  intrinsic  breakdown  mechanism  begins  to  compete  with  avalanche 
ionization,  the  trend  would  be  to  push  the  experimental  points  upward  and 
to  the  left  of  the  prediction  curves  based  on  the  avalanche  effect  alone. 

In  summary,  intrinsic  laser-induced  damage  has  been  shown  to  be 
a  time-dependent  process.  As  the  laser  pulsewidth  was  decreased  to  15  ps 
•he  damage  field  in  NaCl  increased  to  over  107  V/cm.  From  the  pulsewidth 

ferns*  mavRbeS|1mnnohr|Sa  S,u?BC“ted  ‘Private  communicat'ion)  "that "local  field"  of-' 

sgsigliSasgr' 

.  ...  ,  e  ectron  s  rapid  movement  across  the  unit  cell  Tho 

alidity  of  this  latter  assumption  has  not  been  established. 
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dependence  of  the  optical  damage  field,  a  field-dependent  ionization  rate 
was  determined  and  found  to  agree  at  least  qualitatively  with  experiments 
using  dc  fields.  The  agreement  underscores  the  basic  similarity  between 
intrinsic  laser-induced  damage  at  1.06  pm  and  dc  dielectric  breakdown 
and  adds  further  support  to  the  existence  of  avalancne  ionization  in  oelf- 
focusing  "filaments"  ,  where  effects  also  occur  of  a  picosecond  time 
scale. 

We  wish  to  thank  M.  Bass,  L.  Holway,  E.  Yablonovitch  and  J.M. 
McMahon  for  numerous  helpful  discussions.  M.  Bass  and  J.M.  McMahon 
were  particularly  helpful  in  establishing  arrangements  for  this  work. 


37 


V. 


THE  ROLE  OF  INCLUSIONS  AND  LINEAR  ABSORPTION  IN  LASER 
DAMAGE  TO  DIELECTRIC  MIRRORS 


A .  Introduction 

Recent  optical  damage  studies  have  established  intrinsic  limits  for 

9  10  30 

light  propagation  inside  transparent  solids.  *  ’  The  relationship 

between  surface  and  bulk  damage  fields  has  been  determined  in  previous  work 

'4 

by  direct  comparison  of  measured  intrinsic  surface  and  bulk  damage  fields  “ 
and  by  theoretical  analysis.'^  Surface  finishing  techniques  now  exist  which 
for  selected  materials  and  highly  focused  laser  light,  produce  effectively 
imperfection-free  surfaces  whose  damage  fields  are  identical  to  those  of 
the  bulk.-** 

Our  understanding  of  damage  to  optical  coatings,  on  the  other  hand, 
is  not  nearlv  as  advanced  as  our  understandings  of  surface  and  bulk  damage. 
The  reason  tor  this  difference  is  that  material  imperfections  such  as 
structural  defei  ts,  inclusions,  and  residual  absorption  have  not  been 
sufficiently  characterized  in  dielectric  coatings  nor  have  they,  in  fact, 
been  sufficiently  controlled  in  the  manufacturing  process. 

The  present  work  helps  clarify  the  role  of  material  imperfections 
in  dielectric  coating  damage.  By  studying  the  residual  damage  morphology 
and  the  pulse  duration  dependence  of  the  damage  intensity  it  is  shown  that 
damage  thresholds  for  weakly  focused  laser  lights  are  determined  by  highly 
absorbing  inclusions  and  that  there  is  a  large  range  of  inclusion  sizes  pre¬ 
sent  in  the  coatings.  It  is  further  established  that  residual  absorption  in 
coatings  produced  from  transparent  bulk  material  does  not  contribute  to 
optical  damage. 
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1.  Introduction 


By  studying  the  morphologv  of  damage  produced  at  threshold 
and  observing  the  predicted  ^  "pulse  hirst  ion-inclusion  size" 
relationship  we  have  found  that  laser  damage  to  dielectric  coatings  is 
primarily  determined  by  the  presence  of  metallic  or  highly  absorbing 
non-metallic  inclusions.  Specifically,  inclusions  less  than  ~  0.4  pm 
in  diameter  are  most  easily  damaged  by  single  20  psec  duration  ruby 
laser  pulses  while  inclusions  with  diameters  greater  than  ~  3. 5  jini 
determine  damage  resistance  due  to  20  nsec  pulses.  This  suggests 
that  coating  damage  resistance  to  Q-switched  pulses  can  be  improved 
by  eliminating  the  larger  inclusions. 

To  study  linear  absorption  in  films,  experiments  were  performed 
in  which  the  mirror  sample  was  irradiated  by  a  pair  of  pulses  separated 
by  an  interval  of  several  nanoseconds.  Significant  linear  absorption 
would  be  indicated  if  neither  pulse  alone  could  produce  damage,  but 
damage  occurred  when  both  pulses  were  used.  These  experiments 
showed  that  linear  absorption  in  films  does  not  determine  the  damage 
resistance  of  coatings  when  they  are  properly  prepared  from  materials 
which  do  not  show  hulk  absorption. 
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1  xptTi mental  Conditions 


-•  1  Damage  Apparatus 

Experiments  have  been  performed  at  a  wavelength  of  0.69  pm 
using  Caussian-mode  pulses  with  durations  between  20  psec  and  20  nsec. 
The  picosecond  pulses  were  selected  from  a  mode-locked  pulse  train, 
while  the  longer  pulses  were  generated  in  a  dye-Q-switched  oscillator, 
hor  both  systems,  the  temporal  profile,  pulse  energy,  and  a  magnified 
image  of  the  laser  spot  incident  on  the  damage  specimen  are  recorded 
for  each  firing.  Details  of  the  two  systems  are  reported  elsewhere.  ‘5o' 


The  Q-switched  system  has  recently  been  modified  by  the  installa¬ 
tion  of  an  electro-optic  pulse  shaper  between  the  oscillator  and  the  amplifier. 
This  has  allowed  generation  of  strictly  band*,  'th-lirnited  pulses  1.4  nsec 
in  duration  ,  and  pairs  of  pulses  separated  by  a  predetermined  and  re¬ 
producible  interval.  All  damage  experiments  with  shutter-shaped  pulses 
were  monitored  in  the  same  fashion  as  has  been  described  elsewhere.  3  3  *  38 

2.2  Damage  Specimens 

With  a  single  exception  noted  in  the  next  paragraph,  all  damage  speci¬ 
mens  were  electron-gun  deposited,  quarter-wave  stacks  of  either  Ti02/SiC>2 
or  Zr02/Si02  on  fused  silica  substrates.  The  reflectivity  of  each  of  the  mir¬ 
rors  exceeded  90  percent  at  0.69  pm.  A  large  number  of  samples  from  sev¬ 
eral  manufacturers  were  examined;  the  data  presented  is  generally  true  for 
all  specimens. 

A  single  set  of  experiments  was  performed  on  a  ThF^/  ZnS  reflector. 
This  exception  will  be  noted  in  the  text. 


3.  Inclusion  Damage  in  Coatings 


3.  1  Damage  at  20  psec 

Characteristic  "near-threshold"  damage  induced  by  a  20  psec  pulse31**38 
focused  to  a  spot  190  urn  in  diameter  (FWHM  of  the  intensity  distribution) 


'  XhlbltGd  the  Scanning  electron  micrograph  (SEM)  of  Figure  18  * 
T,  e,  ,arSe  ridges  in  the  photo  on  the  left  one  imperfections  in  ^oid 

D  B,aPPl,ed  “  Permil  sca™‘"R  electron  microscopy,  and  are  not 
pt  operttes  of  the  dielectric  coatings  or  the  laser  damage.  The  laser 

::r ,s  thp  array  °f  smm  c  ra,°rs  in  ,he  ^  °f  *«■  paraph 

corresponds  to  the  central  part  of  the  irradiated  region  The 
act  tha,  amage  is  not  produced  in  the  weakly  , eradiated  L  ^  rl g  oes 

“  C°"SlSlent  WUh  th=  «« »<«"■»«",  of  damage  threshold  to  this  ievel  of 
-radiation  (3. 5  J/ cm2  in  ,0  psec.  or  1.7  *  to* *  watts/ cm2, 

the  nhn.A  fied  'icw  °f  “veral  of  the  small  craters  is  shown  in 

o  on  e  right.  Craters  are  randomly  located  in  the  plane  of 
10  coat, „g  but  restricted  primarily  to  the  top  two  layers  of  the  stack. 
There  are  between  2  and  10  x  lo6  such  craters  per  cm2. 

Since  the  size  at  the  top  of  the  craters  may  be  characteristic 
of  explosive  rupture  of  the  material,  a  more  meaningful  measure  of 
the  size  of  the  damage  site  can  be  obtained  at  the  bottom  of  the  crater 

ca™  S‘te  '‘T1"  °btai"ed  USi"B  ‘hiS  CriteHo"  is  <>•*  «».  In  no 
a  e  w„h  a  single  20  psec  pulse  has  a  site  with  diameter  greater  than 
u.  j  u m  heen  observed. 


It  should  be  emphasized  that  the  damage  described  in  this  sec- 
m„  is  caused  by  a  pie  2„  psec  pulsc  and  ,s  nQ,  nccessari,y 

at.c  of  damage  produced  by  a  train  of  several  such  pulses.  The  d„- 
mction  between  these  two  cases  will  be  discussed  in  Section  4. 

3 , 2  Damage  at  1.4  nsec 

The  SEM  in  Figure  10  illustrates  threshold  damage  caused  by 
a  laser  pulse  1.4  nsec  in  duration  focused  to  a  diameter  of  130  pm  at 
e  mirror  sample  surface.  The  energy  density  was  14  J/em2.  Again 

iS  an  array  °f  rand°mI.'  located  small  craters. 

disloi  fed  sinceSth0eWspecfm“n  ?s  "viewed  a t  *n* angl? of  4 TiETd"* 
long  a xis”  ot  J h e'e Hi psl .' ‘ '  meaaurempnts  on  a  diagonal  aligned"™  the 
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onh-  now  the  mean  crater  diameter  is  1 . 5  pm.  The  density  of  these 
i  rs,  10  per  cm  ,  is  less  than  the  density  of  damage  sites  pro¬ 
duced  by  picosecond  pulses.  These  craters  often  extend  furthei  into 
the  coating  than  the  first  one  or  two  layers. 

Smaller  craters,  such  as  those  resulting  from  irradiation  with 
picosecond  pulses,  were  never  produced  by  these  longer  laser  pulses. 

3.3  Damage  at  20  nsec 

The  SEM  in  F  ignre  20  illustrates  threshold  damage  caused  by  a  laser 
pulse  23  nsec  in  duration  focused  to  a  diameter  of  400  pm  at  the  sample  sur¬ 
face.  The  energy  density  was  16  ,T/  cm2.  The  mean  site  diameter  was 

4.3  pm  and  the  site  density  was  104  per  cm2.  In  order  to  measure  so  low 
a  site  density,  it  w,as  necessary  to  irradiate  a  larger  area  of  the  sample. 

t  raters  produced  by  20nsec  pulses  frequently  extend  into  several  layers 
of  the  coating. 

3.4  Interpretation  of  Threshold  Damage 

It  has  been  shown  that  small  impurity  inclusions,  opaque  at  the 
laser  Irequency  and  in  good  thermal  contact  with  the  surrounding  medium 
can  be  heated  by  a  laser  pulse  to  temperatures  in  excess  of  the  melting 
point  of  the  medium.  The  material  near  such  inclusions  can  undergo 

a  phase  change  and  so  laser  damage  results.  When  damage  due  to  the  pres¬ 
ent  e  of  inclusions  occurs,  the  residual  morphology  is  characterized  by  an 
array  of  many  damage  sites  randomly  located  within  the  irradiated  volume.9,  10,  39 
A  similar  morphology  would  be  expected  if  the  damage  were  due  to  any  other 
type  of  localized  material  or  physical  defect.40 

A  morphology  characterized  by  randomly  located  sites  is 
precisely  that  of  the  threshold  coating  damages  shown  in  Figures  18,19, 
and  20.  Vte  conclude  therefore  that  localized  defects  are  the  principal 
cause  of  laser  damage  to  these  optical  coatings. 
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'WHM  in  the  intensity  profile). 


A  hint  that  the  defects  are  impurity  inclusions  is  obtained 
from  the  size  of  the  sites  damaged  by  laser  pulses  of  different  dura¬ 
tions.  By  considering  both  the  heating  and  cooling^  of  spherical,  metallic 
inclusions  in  laser  glasses,  Hopper  and  Uhlmann  and  Bliss  have 
shown  that  for  a  particular  laser  pulse  duration  a  particular  size  of  in- 
elusions  will  be  most  easily  heated  to  damaging  temperatures. 

These  studies  show  that  short  pulses  can  more  easily  damage  small  in¬ 
clusions  than  long  pulses  and  vice  versa.  From  Eq.  4  in  Ref.  11  for  the 
lr  elusion  temperature  as  a  function  of  the  radius,  R^,  one  can  easily  show 
that  to  a  first  approximation  the  inclusion  wnich  is  most  readily  heated  by 
a  laser  pulse  of  duration  r  ,  has  radius  proportional  to  (rp)  .  Thus, 
for  the  three  pulse  durations  used  in  this  work  we  would  expect  the  damaged 
inclusions  to  have  radii  in  the  ratios 


R. (23  nsec):  R  (1.4  nsec):  R.(20  psec)  -  34:8.3:1 

i  i  1 


Note  that  for  a  particular  pulse  duration  a  range  of  sizes  of  "  easily  dam¬ 
aged  inclusions"  can  be  expected  because  the  relationship  between  the  tem¬ 
perature  reached  and  the  inclusion  size  is  not  a  very  sharply  peaked  func¬ 
tion.  11  In  addition,  if  the  level  of  irradiation  is  even  slightly  above  the 
absolute  minimum  required  for  damage  then  several  different  sizes  of  in¬ 
clusions  can  be  heated  to  damaging  temperatures. 


Because  there  is  this  range  of  site  sizes  in  each  case,  the  ratios 
of  the  mean  site  diameters, 

£[.(23  nsec):  Hi<  1.4  nsec):  H.(20  psec)  =  22:7.5:1 

are  to  be  compared  with  the  expected  ratios.  These  sets  of  ratios  are 
in  acceptable  agreement  in  view  of  the  difficulty  in  finding  the  exact  size 
of  the  inclusions  which  produced  the  damage  sites. 

Hopper  and  Uhlmann  have  also  shown  that  for  non- metallic  in¬ 
clusions  having  moderate  absorptivity  there  is  no  "pulse  duration-inclusion 

*  The  treatments  in  Refs.  11  and  34  require  that  the  inclusion  be  completely 
embedded  in  the  surrounding  material.  The  morphology  of  the  damage  in¬ 
dicates  that  this  condition  is  satisfied  for  almost  all  of  the  inclusions  which 
contribute  to  coating  damages. 
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size"  relationship  as  above.  **  However,  if  the  absorptivity  of  the 

inclusion  is  very  high  (i.e.  1  /a  <  R  where  a  is  the  absorption  coeffi- 
_  1  1 

cient  in  cm  at  the  laser  wavelength)  the  same  model  as  used  for  metallic 
inclusions  can  be  applied.  Thus,  the  observed  damage  morphology  and 
"pulse  duration-inclusion  size"  relationship  suggest  that  the  damaging  de¬ 
fects  are  either  metallic  or  very  highly  absorbing  non-metallic  included 
impurities. 

Inclusions  of  either  type  may  arise  from  impurities  in  the  starting 

material,  incomplete  oxidation  of  the  high  index  material  (particularly  in 

41 

the  case  of  TiC^),  introduction  of  material  from  the  electron  gun, 
dust,  or  general  deposition  of  dirt  from  the  chamber.  All  of  these  indicate 
the  possibility  of  improving  the  damage  resistance  of  optical  coatings  by 
improving  the  coating  process  control.  Note  that  long  pulse  damage  re¬ 
sistance  can  be  increased  simply  by  eliminating  the  large  inclusions. 

Non-metallic  absorbing  inclusions  are  likely  to  be  more  strongly 
absorptive  at  shorter  wavelengths  and  in  the  infrared  than  they  might  be 
near  1  Thus,  their  importance  in  damage  problems  will  grow  as  the 

need  for  high  peak  power  devices  in  these  spectral  regions  increases. 

3.5  Above-Threshold  Damage 

The  morphology  of  damage  sites  produced  by  irradiation  at  levels 
significantly  above  threshold  is  consistent  with  the  conclusion  that  damage 
is  caused  by  metallic  or  highly  absorbing  inclusions.  In  Figure  21  there  are 
four  SEM1  s  of  damage  produced  by  single  20  psec  pulses  of  successively 
higher  energy.  These  photographs  are  highly  suggestive  that  damage  for 
above-threshold  irradiation  is  the  result  of  the  compounding  of  damage  at 
many  independent  sites.  Above-threshold  damage  with  single  1.4  nsec  pulses 
appears  to  follow  the  same  pattern. 

A  possible  exception  to  the  observation  that  massive  damage  is 
a  compounding  of  damage  at  many  small  inclusions  occurs  with  the  20  nsec 
duration  pulses.  Above-threshold  damage  sites  at  this  pulse  duration  fre¬ 
quently  consist  of  a  smooth  area,  one  or  two  layers  deep  which  is  similar 
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in  size  and  shape  to  the  incident  laser  beam  (see  Figure  22).  A  number 
of  inclusion  damage  craters  may  occur  either  within  or  outside  the  edges 
of  the  large  site.  Since  the  inclusion  craters  are  seldom  centered  on  the 
larger  area,  it  is  not  clear  that  the  larger  area  results  from  massive  in¬ 
clusion  damage.  In  addition,  the  material  surrounding  these  sites  is  not 
strewn  with  small  inclusion  damages  as  in  Figure  21. 

The  appearance  of  these  large  damage  sites  suggest  that  other  dam¬ 
age  interactions  might  take  place  with  20  nsec  pulses.  Numerous  mechanisms 
involving  combinations  of  inclusion  damage  and/or  intrinsic  processes  can 
be  invoked  to  account  for  the  correlation  between  the  shape  of  the  laser 
beam  and  that  of  the  larger  sites.  One  possibility  is  that  the  outer  layers 
are  simply  heated  to  a  damaging  temperature  due  to  irradiation  of  the  great 
number  of  small  absorbing  inclusions  present.  Linear  absorption  by  the 
coating  materials  will  also  contribute  to  a  damage  morphology  similar  to 
the  laser  beam  cross-section.  This  is  explored  in  detail  in  the  double-pulse 
experiments  described  below. 

It  is  to  be  emphasized  that  observations  concerning  the  potential 
importance  of  other  mechanisms  for  20  nsec  irradiation  in  excess  of  threshold 
in  no  way  effects  the  conclusion  that  threshold  at  this  pulse  duration  is  de¬ 
termined  by  large  inclusions. 

4.  Double-Pulse  Damage  Experiments 

If  a  damaging  quantity  of  energy,  H,is  absorbed  from  a  single 

1.4  nsec  pulse  of  energy, E,  the  same  energy, H, can  be  supplied  by  two 

pulses  each  of  energy  E  (E/2  <  E  <  E)  provided  that  the  interpulse 

P  P 

interval  is  shorter  than  the  cooling  time  of  the  absorbing  volume.  Addi¬ 
tion  of  the  two  pulses  to  produce  damage  can  occur  in  a  uniformly  absorb¬ 
ing  film  or  at  inclusions  too  large  to  be  effectively  heated  by  single  1.4  nsec 
pulses.  This  would  not  occur  at  small  inclusions  which  could  cool  during 
the  interpulse  interval,  or  if  damage  is  due  to  a  fast  response  mechanism 
such  as  an  electron  avalanche  breakdown. 
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Above-Threshold  Damage  in  a  Dielectric  Mirror  Produced  by  a 
23  nsec  Duration  Pulse  Focused  to  a  Spot  Size  of  400/um.  Five 
small  inclusion  craters  are  located  within  the  damaged  area  which 
is  itself  two  layers  deep.  The  area  has  a  shape  like  that  of  the 
incident  laser  beam,  and  is  not  centered  on  the  obvious  inclusion 
damage. 


As  rioted  in  Section  2.  1,  pairs  of  optical  pulses  are  readily 
Renerated  by  the  pulse-shaping  shutter.  The  Pockels  cell  in  this  shut- 
er  is  constructed  as  an  electrical  transmission  line  element.  The 

termin'  Td?6  Pr°Paga*eS  *hr°Ugh  *•“  P°CkdlS  ““  f"‘°  a  cable 
in  a  ref  ,  load.  Removing  the  load  results 

second  ,  °  VOltage  PUlSe’  WhiCh  in  turn  °Pens  the  shutter  a 

nd  time.  Smgle-and  double-pulse  operation  of  the  shutter  are  both 
illustrated  in  Figure  23. 

During  the  experiment,  a  mirror  was  probed  by  single  1.4  nsec 
uses  a,  a  number  of  sites  to  determine  the  energy  E.  required  to  pro- 

duri\.ath',Se  W'th  SmRle  PUlSeS‘  A  given  site  was  irradiated  only  once 
duri.ig  this  sequence.  Additional  sites  were  irradiated  with  a  series  of 

■several  single  subthreshold  pulses  at  intervals  greater  than  one  minute 

■s  ..  general  rule,  even  several  pulses  with  intensities  less  than  that  re: 

quire  or  single  pulse  damage  produced  none.  This  implies  that  perma- 

nen  amage,  undetectable  by  optical  microscopy,  did  not  occur  with  single 

s  Mhreshoid  pulses.  Subsequent  observations  of  damage  due  to  the  addi! 

„  °f  tW°  Cl°SCly  Spaced  could  therefore  be  of  thermal  origin. 

A  record  of  an  experiment  performed  in  this  fashion  is  shown  in 

cafeTh  ,'h  ,  aPPr°Ximate  *hreShoW  for  by  single  pulses  is  indi¬ 

cated  by  the  horizontal  dashed  line,  and  a  lack  of  cumulative  undetected 

damage  by  sequences  of  near-threshold  shots  which  were  fired  on  a  given 
site  without  causing  damage.  6 

A  record  of  some  of  the  double-pulses  experiments  performed  on 

™  TT  ,S  Sh°Wn  ‘n  FigUre  25‘  The  lotal  ““n  (2E  )  of  each  pair 
es  satisfies  the  requirement  2E(  >  2E  >  E  so  that  damage  was 

possible  only  ,f  summation  occurred.  In  several  cases,  damage  did  indeed 
occur  Damage  sites  produced  by  addition  on  this  mirror  were  character- 
is  ic  of  damage  due  to  large  inclusions.  Such  inclusions  would  not  be  ex¬ 
pected  to  damage  upon  irradiation  by  single  1 . 4  nsec  pulses.  This  ob¬ 
servation.  coupied  with  the  failure  of  addition  at  many  sites  is  taken  as 
evi  ence  that  there  is  no  severe  linear  absorption  present  in  this  mirror 
beyond  that  presented  by  the  localized  inclusions. 
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SHAPED  INPUT 

OUTPUT  Fig.  23  Shutter  for  Generating  Single  Pulses  of  1.4  nsec  Duration,  or 

Pairs  of  1.4  nsec  Pulses  Spaced  by  a  Reproducible  Interval. 
Removing  the  load  from  the  cable  causes  a  voltage  reflection 
which  gates  the  Pockels  cell  a  second  time. 


I 
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The  observation  that  large  inclusions  can  be  damaged  due  to 
(.  imulative  heating  by  successive  pulses  can  be  generalized  to  explain 
damage  by  a  train  of  mode-locket  >ulses.  If  the  intensity  of  each  of 
t!  pulses  is  such  that  singly  they  cannot  produce  damage,  damage  to 
large  inclusions  can  still  be  the  principal  failure  mechanism.  This 
mav  explain  the  difference  between  the  damage  threshold  found  for 
single  20  psec  pulses,  ~  10  watts/ cm  ,  and  the  intensity  limits  of 
mode-locked  oscillators  set  by  coating  damage,  ~  1-4  X  10  watts/ cm  . 

It  was  a  general  observation  that  there  was  no  consistent  pattern 
of  pulse-pair  addition  indicative  of  material  absorption  in  the  mirrors 
used  in  this  study. 

In  order  to  ascertain  that  the  double-pulse  experiment  does  in¬ 
deed  detect  absorption,  we  have  used  the  technique  to  study  a  ThF^/  ZnS 
reflector  which  is  known  to  be  linearly  absorptive  at  \  =  0.69  pm. 
Tvpical  double-pulse  data  for  this  mirror  is  shown  in  Figure  26,  where 
the  single-pulse  damage  threshold  is  again  represented  by  the  dashed 
line.  For  this  reflector,  the  expected  consistent  addition  of  pulses  is 
seen  to  occur. 

These  results  indicate  that  if  the  coating  materials  are  non¬ 
absorbing  than  a  combination  of  linear  absorption  and/or  absorption  by 
small  inclusions  cannot  be  invoked  to  explain  results  such  as  those  shown 
in  Figure  22.  It  is  possible  therefore  that  this  type  of  damage  was  caused 
by  an  intrinsic  process.  To  study  this  possibility,  experiments  in  which 
a  spot  size  sufficiently  small  such  that  the  irradiated  area  of  the  coating 
is  not  likely  to  contain  inclusions  which  are  easily  damaged  by  the  pulses 
used  are  in  progress. 

5.  Conclusions 

The  principal  conclusions  of  this  work  are: 
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1.  The  practical  damage  threshold  for  optical  coatings  at 
Aim  is  determined  by  the  presence  of  absorbing  impurity  inclu 


2.  The  predicted  "  pulse  duration-inclusion  size"  relationship 
m  metallic  or  highly  absorbing  non-metallic  inclusions  was  observed. 

3.  \  cry  small  (  <  0.3  Aim)  inclusions  are  the  limiting  factor 
m  coating  damage  resistance  to  single  20  psec  ruby  laser  pulses. 

4.  Larger  inclusions  (  >  4  Aim)  are  responsible  for  threshold 
coating  damage  when  long  duration  pulses  (*  20  nsec)  or  trains  of  mode- 

ed  pulses  are  used.  Damage  resistance  to  such  irradiation  would  be 
improved  if  only  this  class  of  larger  inclusions  could  be  eliminated.  In- 
isions  of  these  sizes  are  appropriate  to  laser  radiation  in  the  visible  or 
n-ar  infrared.  For  other  wavelengths  the  limiting  inclusion  sizes  must 
l,e  scaled  according  to  the  wavelength. 


5*  Linear  absorPti°n  in  films  does  not  determine  the  damage 
res.stance  of  coatings  when  they  are  properly  prepared  from  materials 
winch  do  not  show  bulk  absorption. 
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VI. 


RELATIONSHIP  BETWEEN  SELF -FOCUSING  AND  OPTICAL 
DAMAGE 


A .  Introduction 

A  number  of  researchers  have  analyzed  the  self-focusing  effect 
by  using  a  paraxial  ray  formula.  This  approach  has  been  found  to  be 
inadequate  in  the  ultra-high-power  regime  where  small-scale  beam 
distortions  can  lead  to  effects  that  are  absent  at  lower  laser  powers. 
However,  the  paraxial  ray  model  has  been  useful  in  describing  both 
incipient  self-focusing  at  moderate  laser  powers  and  the  competition 
between  self-focusing  and  diffraction  at  low  laser  powers  near  the 
critical  power  P(>^  for  catastrophic  beam  collapse. 

The  relationship  between  self-focusing  and  optical  damage  is 
discussed  in  this  section  within  the  framework  of  a  paraxial  ray  model. 
Because  of  the  low  laser  powers  considered  (P  <  P^),  the  paraxial 
ray  model  appears  to  provide  an  adequate  description  of  beam  propaga¬ 
tion.  The  analysis  summarizes  the  effects  of  the  index  nonlinearity  on 
a  focused,  low-power  laser  beam  and  shows  how  this  effect  can  be  used 
to  measure  self-focusing  parameters.  As  an  illustration  of  the  appli¬ 
cation  of  this  measurement  technique,  the  self-focusing  parameters  for 
sapphire  are  determined. 
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Abstract 

A  new  method  for  measuring  the  self-focusing  parameters  of 
transparent  materials  is  reported.  The  technique  consists  of  using 
lenses  with  different  focal  lengths  to  induce  intrinsic  bulk  damage  at 
low  powers  where  catastrophic  self-focusing  is  absent  and  then  com¬ 
paring  the  results  to  analytical  calculations.  Using  this  technique,  the 
self-focusing  parameters  for  sapphire  are  determined. 
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It  has  been  observed  that  a  third-order  optical  nonlinearity  can 

lead  to  catastrophic  self-focusing  of  light.42  This  nonlinearity  causes 

the  index  of  refraction  to  increase  by  an  amount  n.,<E  /  in  the  presence 

of  an  optical  electric  field.  Values  of  n2  for  various  materials  can  be 

determ, ned  from  measurements  of  induced  birefringence,  43  by  studies 

of  the  movement  or  positions  of  optical  damage  tracks  in  crystals, 

or  by  observing  the  changes  in  the  intensity  distribution  ol  a  transmitted 

laser  beam  as  the  input  power  is  increased.46  The  electric  contribution 

to  n2  can  also  be  inferred  from  measurements  of  optical  third  harmonic 
generation. 


Reported  here  is  a  new  technique  for  measuring  self- focusing 
parameters  consisting  of  determining  the  powers  necessary  to  induce 
intrinsic  bulk  damage  as  lenses  with  different  focal  lengths  are  used  to 
focus  the  radiation  and  compare  the  results  to  analytical  calculations. 
High  input  powers  are  not  needed  so  that  relatively  small,  stable  lasers 
can  be  used.  This  technique  is  described  below,  and  as  an  illustration 

of  its  application,  the  self-focusing  parameters  for  sapphire  are  deter¬ 
mined. 


We  consider  a  low  power  laser  beam  with  input  power  P  <  P 
where  P^  in  esu  units  is48-50  cr 

Pcr  =  K02  c/32  n2  (1) 

The  vacuum  laser  wavelength  is  and  c  is  the  vacuum  light  velocity. 
For  such  powers  catastrophic  self- focusing  is  absent,  and  the  focused 
beam  diameter  is  decreased  from  its  diffraction  limit  of  dQ  to  a  smaller 
value,  dj,  because  of  the  index  nonlinearity.  As  a  result  of  this  reduc¬ 
tion  in  beam  diameter,  the  on-axis  intensity  is  increased  by  the  factor 

(d0/dl)  ’  a  factor  which  is  a  function  of  both  the  laser  input  power  and 
the  nonlinear  index  n2»  (See  Fig.  27) 
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-  Light  ray  paths  with 

Light  ray  paths  with 

and  R  <  P 

input  cr 


n2  *0 
n2  #0 


Fig.  27 


The  Effect  of  the  Index  Nonlinearity  on  the  Diameter  of  the  Focal 

Wh  "  «  ,'”Cy  10  self-focus  acts  in  opposition  to  diffraction. 

Whereas  diffraction  produces  l  positive  ray  curvature,  self- 
focusing  produces  a  negative  ray  curvature.  The  focal  diameter 
inside  a  nonlinear  medium  (n,  *  0)  will  be  effectively  reduced, 
therefore,  from  dQ  --  its  value  with  n2  ■  0  --  to  di,  and  the  on- 
xis  light  intensity  at  the  focus  will  be  correspondingly  increased 
s  shown  by  Wagner  et  al.49,  the  presence  of  the  index  nonlinearity 
also  causes  the  point  of  minimum  beam  radius  to  be  displayed  to 
the  right  toward  the  geometrical  focus.  P  y  ° 
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lculations  jf  d0/dj  have  been  made  for  Gaussian  beams  from  a 
oa  y-state,  paraxial  ray  analysis.  These  calculations  show  that  if  r 
ts  the  on-axts  tntensity  of  a  strongly  focused  laser  beam  when  =  0 
then  the  on-axis  intensity  I,  in  the  presence  of  the  index  nonlinearity  is* 

'd  =  !0  (I  -  p/pcr>  •  (2) 

As  the  input  power  is  increased  to  within  a  few  percen.of  P  corrections 
Lq.  (2)  may  become  important48  and  numerical  calculations  may  be 
required  to  accurately  establish  the  relationship  between  Id  and  IQ. 

The  steady-state  analysis  summarized  by  Eq.  (2)  is  appropriate 
or  tightly  focused  laser  beams,  because  the  electronic  contribution  to 
its  steady-state  whenever  the  laser  pulsewidth  t  is  greater  than  a 
response  time  r  given  by  the  ratio  of  the  -radius  of  theP focused  laser  beam 
and  the  radial  compression  wave  velocity  v  51  I„  sapphire  r  =  lo'8  sec 
when  the  focused  beam  diameter  is  220  pm  and  the  propagation  direction 
r  along  the  crystalline  c-axis.  If  elec, restriction  is  slightly  transient. 

Eq.  (1)  may  be  used  with  Pc_.  corrected  in  a  manner  discussed  by  Kerr.52 

Having  this  determined  the  intensification  produced  by  the  index 
nonlinearity,  we  look  for  a  process  which  is  sensitive  to  a  constant 
value  of  on-axis  laser  intensity.  The  powers  necessary  to  reach  this 
in  ensity  will  be  recorded  as  dp  (see  Fig.  27)  is  varied  by  using  lenses 

wi  .  erent  focal  lengths.  Intrinsic  optical  damage  from  avalanche 

reakdown  is  a  convenient  choice  for  this  intensity-sensitive  process.9'10 
though  laser-induced  avalanche  breakdown  is  an  inherently  statistical 
process,  ,t  is  virtually  threshold-like  for  many  materials.  A 

damage  intensity  Id  can  be  defined10  as  tha,  value  of  intensity  inside  the 
material  which  is  necessary  in  order  to  induce  damage  on  a  single  shot 
with  a  probability  of  0.  5.  The  damage  power  is  the  incident  power 
necessary  to  reach  Id  under  particular  external  focusing  conditions. 
Breakdown  leads  to  a  spark  and  local  melting  of  the  material  and  is 
therefore  easily  observed. 

Prof.  J.  Marburger  has  pointed  out  thif  F'r,  ~  . - 

result.  Using  Eqs.  (2  If!)  anri  <7  approximate 

bee^focu^d^iTfficiently^to'induce'kitrinsi'^opticarda'mag^T 


Pcr  is  determined  from  measurements  of  the  damage  power  P 
made  with  various  lenses  having  different  focal  lengths.  Complications 
from  aberrations  in  the  focusing  optics52  and  inaccuracies  in  deter¬ 
mining  the  input  powers  can  be  greatly  reduced  by  the  use  of  a  reference 
crystal  such  as  NaCl  which  can  be  damaged  over  a  large  range  of  input 
powers  before  self- focusing  becomes  important.10  Approximate  values 

of  Pcr’  which  are  useful  for  the  initial  selection  of  focusing  optics,  can 
be  calculated  theoretically.  47,  ,52 

To  analyze  the  data,  it  is  convenient  to  rewrite  Eq.  (2)  as 

l/P  =  (2n0/H  !-lj  +  1/P  (3) 

w 

Here  we  have  assumed  that  the  intensity  distribution  at  the  focus  in  the 
absence  of  the  index  nonlinearity  is  I(r)  =  IQ  exp  (-2r2/w2)  where  w  =  d  /2, 
and  we  have  replaced  IQ  in  Eq.  (2)  by  2nQP/rw2.  The  focal  lengths  o?  thl 
various  lenses  can  be  related  to  the  corresponding  values  of  w.  The 
damage  intensity  Id  does  not  depend  on  w  for  avalanche  breakdown.  10 
A  plot  of  1/P  as  a  function  of  l/w2  then  yields  a  straight  line  with  l/P 
as  the  intercept  and  2nQ/ld  as  the  slope.  Cr 

Two  precautions  should  be  taken  when  using  this  measurement 
technique.  The  first  is  to  avoid  damage  from  absorbing  inclusions,  and 
the  second  is  to  account  for  any  dependence  of  avalanche  breakdown  on 
the  laser  pulse  width.  Methods  for  distinguishing  between  inclusion  and 
intrinsic  damage  are  detailed  in  Refs.  10  and  18.  The  laser  pulsewidth 
dependence  to  the  damage  field  has  been  measured.17’54  * 


As  an  illustration  of  the  use  of  this  technique,  we  consider  data 


accumulated  in  two  different  experiments  on  optical  damage  in  sapphire. 
In  the  first  experiment53  a  Q-switched  TEMQ0  ruby  laser  was  used  by 


Giuliano  et  al.  to  damage  sapphire  with  three  different  lenses.  In  the 

Li  x 


second 


a  measurements  of  Id  was  made  with  two  different  lenses  and 


An  increase  in  damage  intensity  of  no  more  than  about  15  percent  was 
predicted  as  the  laser  pulsewidth  is  varied  from  20  to  5  ns.  This 
prediction  was  confirmed  in  Refs.  51  and  24.  A  frequency  dependence 
to  the  breakdown  intensity  will,  if  present,  act  to  compensate  this 
° Pecause  of  Possible  inaccuracies  introduced  from  the  pulse¬ 
width  and  frequency  dependence  to  breakdown,  experiments  designed 

£  th°  teC^TuqUe  r1eport?d  in  the  text  should  be  conducted  at  a  single 
frequency  with  the  pulsewidth  being  held  reasonably  constant. 
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a  Q-switched  TEM()0  YAG:Nd  laser.  The  pulse  widths  were  20  nsec  and 

5  nsec,  respectively.  54  Table  III  summarizes  the  relevant  results  of 
these  studies. 


In  order  to  compare  the  data,  the  wavelength  dependence  of  P 

must  be  considered.  It  can  be  seen  from  Eqs.  (2)  and  (3)  that  the  data** 

can  be  plotted  together  if  the  coordinate  axes  for  the  YAG  data  are 

properly  scaled  relative  to  those  for  ruby.  In  the  plot  of  l/P  versus 

1/w  in  Fig.  28,  the  axes  are  scaled  by  f\Q  (YAG)/\  (ruby)  ]  2  =  2.35. 

The  intercept  on  the  vertical  ruby  axis  is  Pcr  at  0.6943  pm,  and  the 

slope  measured  relative  to  either  pair  of  axes  is  (2/71!  )n 

'  d  0* 

Pcr  1S  found  t0  be  °*83  X  !06  watts  at  ruby  wavelength;  n9  = 

0.56  X  10-13  esu;  and  Id  is  1.1  X  1011  watts/cm2.  The  value  of  n9 
hes  within  the  range  of  values  reported  in  Ref.  44  and  I  agrees  with 
the  results  of  Ref.  24. 

2  The  PlQt  in  Fig*  28  curves  downward  when  1/w2  £  0.02  X  10-2 
(pm)  because  the  electrostrictive  contribution  to  the  critical  power  is 
becoming  transient.  Even  though  Eq.  (3)  is  not  valid  for  transient  ef¬ 
fects,  we  can  modify  it  to  approximate  slightly  transient  behavior  by 
letting  Pcr  m  Eq.  (1)  become  a  function  of  w  as  described  by  Kerr.52 
If  it  is  assumed  that  the  electrostrictive  nonlinearity  is  much  larger  than 
the  electronic  effect,  a  damage  power  of  about  1.  1  X  106  watts  is  there¬ 
by  predicted  for  w  =  140  pm.  This  value  compares  favorably  to  the  meas¬ 
ured  value  of  1.2  X  10  watts  in  Table  III. 

The  electronic  contribution  to  n2  can  in  principle  be  determined 
by  the  present  technique  if  electros triction  is  negligible  because  conditions 
of  measurement  have  been  chosen  so  that  t  h  <<  1.  Because  of  the  re¬ 
quirement  that  the  input  power  be  less  than  P^,  t  /r  <<  1  cannot  be 
achieved  by  increasing  t  even  though  the  use  of  weakly  focusing  optics. 

This  extreme  transient  limit  for  electrostriction  can  be  achieved,  how¬ 
ever,  by  using  single,  mode-locked  laser  pulses  with  time  durations 
of  the  order  of  10  sec.  Such  pulses  must  be  temporarily  smooth  and 

hence  transform-limited  so  that  a  YAG:Nd  system  is  a  natural  choice 
for  the  laser. 
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TABLE  III 


DAMAGE  POWER  IN  SAPPHIRE  FOR  DIFFERENT  BEAM  SIZES 

Wavelength  »  ^ml  P  — — 

0.6943  nm  20  0.51  ±.04 

(ref.  15)  53  0.74  ±.07 

140  1 . 23  ±  . 10 

1.064  (im  0.24  ±.05 

(ref.  16)  “•>  0.05  ±.02 


both  experiments  Gaussian  laser  beams  were  used. 
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I/P  (10  6  welts' 


PBN-74-2 14 


1/w2  [10‘2  (Mm)"2] 
2  4 


®  Data  taken  at  1.06 /i,m;  use 
scales  above  and  to  the  right. 

+  Data  taken  at  0.69/nm;  use 
scales  below  and  to  the  left. 


.15  .20 


Fig.  28 


1/w2  [lO'2  (Mm)*2] 

Jhe  Dependence  of  Damage  Power  for  Sapphire  on  Beam 
Focal  Diameter.  A  plot  of  l/P  as  a  function  of  l/w2  produces 
a  straight  line  according  to  Eq.  (3).  P  is  the  input  damage 
^nd  w  is  the  Gaussian  beam  radius.  For  large 
w(l/w  ~  0.02  X  10-2  [|im]  2),  departures  from  a  linear 
relationship  are  expected  because  electrostrirtion  becomes 
transient.  This  is  experimentally  observed  as  seen  in  the 
insert  where  the  scales  for  the  ruby  data  have  been  expanded. 


.SUOM  01)  d/I 


It  should  be  noted  that  the  analysis  summarized  by  Eqs.  (2)  and  ( 3 ) 

is  based  on  the  assumption  of  an  infinite  self-focusing  medium.  In  fact, 
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the  tendency  to  self-focus  begins  at  the  sample  entrance  face.  A  correc¬ 
tion  to  Eq.  (2)  that  accounts  for  the  finite  medium  has  been  evaluated  using 
the  method  of  Ref.  56.  v  It  was  found  that  the  error  introduced  by  the 
assumption  of  an  infinite  medium  is  less  than  10  percent  provided  that  the 
distance  Az  between  the  beam  focus  and  the  entrance  face  satisfies 

1 

Az  >  5.4  X  10'3  (f2>  .  (4) 

cr 

The  lens  focal  length  f  and  Az  are  expressed  in  millimeters.  Equation  (4) 
appears  to  have  been  easily  satisfied  in  the  experiments  of  Refs.  53  and 
24. 


In  summary,  a  new  technique  for  measuring  self-focusing  para¬ 
meters  in  transparent  solids  has  been  developed  which  uses  optical  damage 
induced  by  low-power  lasers.  This  technique  has  been  applied  to  previously 
reported  measurements  of  breakdown  intensities  in  sapphire  to  infer  a  value 
for  n,;. 


The  author  wishes  to  acknowledge  numerous  helpful  discussion 
with  M.  Hass ,  J.H.  Marburger  and  Y.R.  Shen. 


*  Equation  (2-9)  in  Ref.  56  was  integrated  from  zq  to  z  instead  of  -  oo  to 
z.  Here  zQ  is  the  coordinate  of  the  sample  entrance  surface. 
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VH.  INCLUSION  DAMAGE  IN  ZnSe 


r( 


A.  Introduction 


This  section  summarizes  an  analysis  of  inclusion  damage  and  a 
short  experimental  study  of  10.6  pm  laser  damage  in  ZnSe  grown  by 
chemical-vapor-deposition.  The  analysis  is  based  on  a  model  of  in¬ 
clusion  damage  in  which  the  probability  of  finding  a  damaging  inclusion 
inside  the  high  intensity  portion  of  the  laser  focal  volume  is  calculated. 

A  statistical  aspect  to  inclusion  damage  is  predicted,  and  an  operational 
threshold  (probability  i  point)  is  calculated  in  terms  of  the  damage  in¬ 
tensity  of  an  individual  inclusion  and  appropriate  beam  parameters. 

Data  is  presented  that  describes  10.6  pm  laser  damage  in  ZnSe  induced 
by  external  focusing  with  lenses  having  two  different  focal  lengths.  The 
damage  intensity  was  lower  with  the  longer  focal  lens  as  expected  from 

an  inclusion  damage  mechanism.  This  observation  is  quantitatively 
explained. 


B.  LASER  DAMAGE  FROM  INCLUSION  ABSORPTION  IN  ZnSe* 


David  W.  Fradin  and  Dominic  P.  Bua 
Raytheon  Research  Division 
Waltham,  Massachusetts  02154 


ABSTRACT 


The  inclusion  damage  studies,  conducted  at  10.6  um  in  ZnSe, 
show  a  dependence  of  bulk  inclusion  damage  levels  on  the  degree  of 
external  beam  focusing.  This  same  effect,  previously  seen  and  ex¬ 
plained  in  thin  film  damage  studies,  is  quantified  for  bulk  inclusion 
damage  by  using  a  simple  model  based  on  a  statistical  approach  to 
inclusion  damage. 
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Hartford,  Conn.  06108. 
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1  •  Introduction 


When  optical  damage  results  from  inclusion  absorption,  the  damage 
intensity  will  depend  on  the  focused  beam  diameter.  60  This  dependence 
exists  because  the  inclusions  responsible  for  damage  are  randomly  dis¬ 
persed  in  the  material.  When  the  beam  diameter  2uQ  is  of  the  same  order 
of  magnitude  as  the  average  inclusion  separation  dQ,  this  dependence  is 
pronounced  and  can  be  calculated  from  the  probability  of  illuminating  a 
damaging  inclusion  with  the  high  intensity  portion  of  the  light  beam.60, 01 
In  the  limit  of  sufficiently  large  beam  diameters  (2uQ  >>dQ)  or  sufficiently 
small  diameters  (2uQ  <<  dQ),  however,  the  damage  intensity  will  be  in¬ 
dependent  of  beam  size.  Constant  damage  intensity  results,  in  the  former 
case,  from  the  fact  that  the  probability  of  intercepting  an  inclusion  in  the 
intense  portion  of  the  light  beam  approaches  unity  for  large  beam  diameters. 
In  the  latter  case,  the  damage  intensity  becomes  independent  of  Up  because 
inclusions  will  be  missed  by  most  tight  focused  laser  pulses  and  damage 
will  thus  normally  result  from  an  intrinsic  process. 

In  Sec.  2  analytical  expressions  will  be  derived  for  the  probability 
of  bulk  inclusion  damage  induced  by  Gaussian  laser  beams.  These  ex¬ 
pressions  will  allow  us  to  quantify  the  ideas  just  expressed  and  to  inter¬ 
pret  the  10.fi  M-m  ZnSe  damage  data  presented  in  Sec.  4. 

2.  Theory 

In  order  to  more  accurately  clarify  the  inclusion  damage  problem 
considered  here,  the  following  questions  are  posed:  If  optical  damage 
from  inclusion  absorption  develops  when  the  laser  beam  is  focused  to  an 
ox-axis  intensity  Iq  ,  what  is  the  threshold  intensity  1^  necessary  to 
damage  an  individual  inclusion?  How  does  1^  vary  with  the  degree  of 
external  focusing? 

Two  simplifying  assumptions  will  be  made  in  the  analysis.  First, 

IQ  will  he  assumed  to  be  less  than  the  intrinsic  breakdown  intensity  I,  . 
Second,  all  the  inclusions  are  assumed  to  be  identical  in  the  sense  that 
each  has  the  same  damage  intensity  I,  when  the  laser  pulse  duration  is 
held  constant. 
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This  second  assumption,  although  it  might  appear  at  first  to  he 
an  unphysieal  restriction,  can  be  a  useful  approximation  to  the  inclusion 
damage  problem.  As  shown  by  Hopper  and  Ulilman,  ^  the  damage  in¬ 
tensity  for  highly  absorbing  inclusions  is  lowest  when  the  thermal  dif¬ 
fusion  time  t  and  the  pulse  duration  are  comparable.  Since  t~R^/I) 
where  R  is  the  inclusion  radius  and  D  is  the  host  thermal  diffusivity, 
threshold  damage  is  principally  determined  by  inclusions  having 


R 


(Dt  )■ 
P 


(1) 


These  inclusion  have  a  well-defined  damage  intensity  Ij.  At  longer  pulse 
durations,  larger  inclusions,  if  present,  will  be  responsible  for  threshold 
damage  and  Ij  will  scale  as  (tp)”2  as  long  as  (1)  can  be  satisfied. 

We  will  see  that  even  when  1^  is  a  well-defined  quantity,  inclusion 
damage  can  appear  to  be  a  statistical  process  because  of  the  random  place¬ 
ment  of  inclusions  inside  the  solid.  For  this  reason  it  is  necessary  to  de¬ 
fine  an  operational  damage  threshold  as  that  intensity  at  which  the 
probability  for  damage  is  This  same  convention  was  used  in  the  ZnSe 
experiment  in  Sec.  3.3  as  well  as  in  intrinsic  damage  studies  where 

1 8  20 

statistics  result  from  the  physical  nature  of  the  damage  process  itself  ’ 
and  not  from  the  inhomogeneous  nature  of  the  medium  as  in  the  present 
case. 


To  derive  the  bulk  damage  probability,  we  begin  with  the  two- 

dimensional  problem  of  points  randomly  distributed  on  a  plane. 

60 

DeShazer,  et  al.  have  considered  this  same  starting  point  in  their 
approximate  calculations  of  the  damage  intensity  of  thin  film  optics. 


Analysis  shows  that  the  probability  that  no  points  will  be  within 
a  circle  or  radius  r  on  the  plane  is^ 


P0  = 


2 

-o  7T  r 

e 


(2) 


where  o  is  the  average  surface  density  of  points.  Consider  a  Gaussian 
light  beam  that  is  incident  on  this  plane  and  that  has  the  intensity  distri¬ 
bution 
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I(r,  z) 


exp 


2r2 

u(z)2 


(3) 


where  w  (z)  [  1  -  (z/zR)“  J  2  and  zR  is  the  Rayleigh  distance  7ruQ2/\  . 

If  each  of  the  points  of  the  plane  corresponds  to  an  inclusion  with  damage 
intensity  Ij,  then  the  probability  that  optical  damage  will  occur  on  this 
plane  is  just  the  probability  that  an  inclusion  will  lie  with  the  region 
where  I(r,z)  >  Ij.  This  region  is  bounded  by  a  circle  of  radius  r^  de¬ 
fined  by  the  intersection  of  the  surface 


J(rT,z)  =  I 


(4) 


with  the  plane  at  constant  z. 


From  Eqs.  (3)  and  (4),  r,p  is  found  to  be 
rT  =  W0  V  in  +  (z/2R)Z]jfn  7/[  1  +  (z/zR)^]  (5) 

whe.-e  y  =  Iq/Ij  >  1.  The  probability  Pp  that  damage  will  occur  is  just 

PD  (I0,W0)  =  1  "  Po  =  1  "  exP  (-7rcrT2)  .  (6) 

If  the  plane  is  located  at  the  beam  focus  (z  =  0),  then  rT  simpli¬ 
fies  and  Pq  can  be  expressed  as 


PD(I0,UJ0) 


(7) 


Here  n  lsjhe  average  number  of  inclusions  inside  the  circle  of  area 
*uq  /2.  n  and  are  related  by 


Equations  (7)  and  (8)  give  the  probability  that  a  focused  Gaussian  beam 
will  induce  damage  on  the  surface  of  a  solid  or  in  a  thin  film  layer. 
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The  analysis  can  row  be  extended  to  three  dimensions.  Before 
calculating  PQ  exactly,  consider  an  idealized  gaussian  beam  that  re¬ 
mains  perfectly  collimated  over  a  distance  For  this  geometry 

PD  is  still  given  by  Eq.  (6),  but  a  is  replaced  by  pjfeff.  The  volume 
density  of  inclusions  p  is  treated  as  an  exogeneous  variable  that  has 
been  determined  by  an  independent  technique  such  as  the  light  scatter¬ 
ing  measurement  described  in  the  next  section-  With  this  substitution 

PD  =  1  ^effrT2)  (9) 

=  1  -  exp  (-p  Veff)  . 

Veff  is  that  portion  of  the  focal  volume  in  which  the  laser  in¬ 
tensity  is  high  enough  to  cause  inclusion  damage. 

For  a  true  focused  beam,  P^  still  has  the  form  of  Eq.  (9).  The 
basic  form  of  the  probability  is  unchanged  for  complicated  geometries 
because  of  the  exponential  behavior  of  pQ  and  the  multiplicative  nature 
of  compound  probabilities.  V  ff  can  be  approximated  by  cones  of  re¬ 
volution  as  noted  in  Fig.  29.  The  radius  of  the  base  of  either  cone  is 
r rj-,  and  the  cone  height  is  z^,  defined  by 

1(0,  zT)  =  Ij  (10) 


or 

zT  =  ZR  (7  -1)2  =  -7^SL-  (7  -  1>*.  (11) 

V  ..becomes,  in  this  approximation, 
ett 

Veff  (cones)  =  n  r^,“  (Q)  Zj,  (12) 

where  r ^,(  0)  is  given  by  Eq.  (5)  with  z  =  0  .  Alternately,  can  be 

calculated  exactly  as  we  now  show. 
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Fig.  29  Approximate  Focal  Volume  for  a  Gaussian  Beam. 
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veff  is  the  cylindricallv  symmetric  region  hounded  by  the  surface 

ol  Eq.  (4).  Making  use  of  the  cylindrical  symmetry,  V  can  be  written 
as 


V 


eff 


■  *1 


zT 

0 


7r  r^,“  (z)  dz. 


(13) 


Inserting  Eqs.  (5)  and  (11)  into  V  ff  gives 

ZF 


V 


eff 


2  (  F  ,  y 

-  J  [  1  +  (z/zR)  ]j?n 


l/[  1  +  (z/zR)  ) 


J  dz. 


(14) 


This  integral  can  be  expressed  in  closed  form  using  standard  integral 
tables  with  the  result  that 

Veff  =  (7ru0  *  ZF  ^  2/9  ^  *7  +  5>  "  tan"1  (7  -  l)2  ]  (15) 

Since  pQ  =  <-pVeff)  and  PQ  =  1  -  pQ  from  Eq.  (6),  it  follows  that 


D 


1  -  exp  (-pVeff) 


1  -  .exp  {-p 


i  2,2 

(7TCJ  q  ) 


^  f-g  (7- D2  (t+5)  - tan'1  (T- I)2  ]}  .  (16) 
As  before,  7  -  Iq/I^,  and  only  the  region  7  >  1  is  considered. 

Equation  (16)  expresses  the  probability  that  inclusion  damage  will 
occur  in  a  solid  and  shows  how  this  probability  depends  on  focused  beam 
radius  (uq),  the  on-axis  intensity  at  the  focus  (Iq),  the  damage  intensity 
for  an  individual  inclusion  (I  ^ ) ,  and  the  measured  inclusion  density 

(p  -  -Aj03>. 

pq  is  not  a  step  function  except  in  the  limit  cj  -  00.  There 
exists  a  range  of  conditions  in  which  0  <  PD  <  1  so  that  optical  damage 
will  appear  to  be  a  statistical  process  as  different  damage  sites  are 
probed.  I  nlike  the  statistics  resulting  from  avalanche  breakdown1^  ’ 2<1 
however,  PD  is  a  function  of  laser  intensity,  not  peak  field,  and  the 
damage  intensity  is  well-defined  at  each  damage  site.  Statistics  are 
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apparent  in  the  inclusion  problem  because  each  damage  site  is  different. 
For  a  perfect  single-mode  laser  held  to  constant  peak  intensity,  each 
site  will  damage  on  the  first  shot  or  no  at  all.  Referring  to  Eqs.  (7) 
and  (9),  it  is  seen  that  the  intensity  range  over  which  PD  is  measurably 
different  from  0  or  1  depends  on  the  average  number  of  inclusions  n 
inside  the  focal  volume.  When  n  is  large,  as  expected  for  weakly  focused 
laser  beams,  PD  approaches  a  step  function  with  the  discontinuity  located 
at  IQ  =  Ij.  The  dependence  of  PD  on  n  is  illustrated  schematically  in 
Fig.  30. 

Because  the  statistical  nature  of  inclusion  damage  precludes  the 
observation  of  a  well-defined  threshold  for  focused  beams,  an  effective 
damage  threshold  Iq^  must  be  defined.  ^  We  let  Iq^  be  that  value  of  Iq 
which  induces  damage  with  a  probability  of  When  Pn  =  |  is  inserted 

D  u 

into  Eq.  (16),  an  equation  results  which  relates  Iq  to  1^.  From  Eq.  (13) 
we  find 

i  i  1  q\ 

(y  -  l)2  (7+  5)  -  6  tan  (7  -  l)2  =  - TT~rl  (10 

2  p(  7T  (J  q  ) 


The  right  hand  side  of  this  equation  can  be  evaluated  in  terms  of  known 
parameters  giving  an  equation  for  7  =  Iq^/I^. 

Equation  (17)  is  our  basic  result.  If  cJq,  p  and  1^  are  measured, 
then  Eq.  (17)  allows  us  to  calculate  1 1 ,  the  damage  intensity  for  a  single 
inclusion.  Having  obtained  Ij,  it  is  then  possible  to  predict  the  damage 
intensity  measured  under  different  conditions  of  focusing  using  Eq.  (17) 
and 


_  2  X  . 
w  0  "  T  D  f* 


(18) 


9 

D  is  the  gaussian  l/e  diameter  at  the  lens  of  focal  length  f.  Alternately, 


D 


p  can  be  calculated  from  two  or  more  measurements  of  Iq  made  with 
different  lenses. 


I 

h 


the  Relationship  Between  P 
erent  Values  of  "n. 


When  7  approaches  unity  (weak  focusing),  Eq.  (17)  can  be  approxi¬ 
mated  by 

u,  _  i>3/2  _  9  £n  2\ 

{y  l)  ^  (19) 

4p(7TU0  ) 

In  the  next  section  we  will  use  Eq.  (17)  to  analyze  the  results  of  inclusion 
damage  studies  in  ZnSe. 

3.  Intrinsic  damage  statistics 


^'llei8ntrinS^C  ^ama^e  Process  has  been  found  to  have  a  statistical 
character.  Although  these  statistics  are  less  pronounced  than  the 
statistics  associated  with  inclusion  damage,  they  should,  in  principle, 
be  included  in  a  complete  model  of  damage  statistics. 

To  illustrate  the  manner  in  which  intrinsic  and  inclusion  statistics 
are  combined,  we  calculate  the  total  damage  probability  P  from  the 
inclusion  damage  probability  Pj^C  and  the  intrinsic  damage  probability 
PD  r*  PD  is  Siven  by 

“  1  "  (probability  that  no  damage  occurs) 

-  i  -  a  -  pj?0)  (i  -  p^1-)  .  (20) 

We  have  used  the  fact  that  the  probability  that  no  damage  occurs  is  given 
by  the  compound  probability  that  damage  does  not  occur  in  either  an  in¬ 
trinsic  or  an  inclusion  process.  Expanding  Eq.  (20)  gives 

P  =  fpinc  4  pintr.  inc  ^.intr 

Pti  PD  4  PD  >  “  PD  PD  *  (2D 

In  the  limit  where  P^  tr  is  a  step  function  defined  as 

In-  T  <  Tintr 

pintr  _  ,,-intr  >  _  J  u*  x0  1 

D  'lo>  '  U. „*.*»*  ’  (22> 
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Equation  (21)  reduces  to 


ptoc.  Io  <  ,in.r 

“  "  U:  10  •  Iin,r  ‘  1231 

Since  it  has  been  found  experimentally  that  Eq.  (22)  is  nearly  valid 
in  ZnSe  for  the  conditions  of  ou-  measurement,  Eq.  (23)  will  be  used  for 

Ptf  Wlth  PDC  given  ^  E(^‘  (16)*  In  experiments  where  the  intrinsic  limits 
are  approached,  however,  it  is  necessary  to  explicitly  consider  the  intrinsic 

damage  probability  in  order  to  adequately  describe  the  statistical  aspects 
of  laser  damage. 


Experimental  measurements  of  inclusion  damage 


ZnSe,  grown  by  a  chemical  vapor  deposition  process  at  Raytheon 
Company,  was  probed  with  laser  light  at  1.06  pm  (t  =  10  ns)  and  10.6  pm 
(tp  =  140  ns)  from  Q-switched  lasers.  The  TEMQ0  outputs  from  both 
these  lasers  were  examined  with  pinholes,  as  described  in  Ref.  10  and 
found  to  be  gaussian  in  form  to  r  ^  3uQ. 

Damage  studies  at  1.06  pm  were  conducted  to  establish  intrinsic 
breakdown  limits  and  to  thereby  provide  a  reference  for  additional  damage 
studies.  Because  the  1.06  pm  beam  could  be  highly  focused  (2W()=  16  pm), 
it  was  possible  to  avoid  inclusions  at  most  damage  sites.  When  inclusion 
damage  did  occur  at  individual  damage  sites,  it  was  possible  to  distinguish 
between  inclusion  and  intrinsic  damage  by  studying  the  residual  damage 
morphology.  ’  Data  corresponding  to  inclusion  damage  were  discarded 
in  the  measurements  at  1.06  pm.  Catastrophic  self-focusing  was  precluded 
by  t  ie  experimental  design9,  10  and  tests  were  conducted  to  confirm  the 
absence  of  measurable  beam  distortion  from  either  a  steady-state  or  a 
transient  index  nonlinearity. 10  In  this  manner  the  incident  intrinsic 
damage  intensity,  I2,  was  measured  at  1.06  pm  to  be  1.2  ±  0.3  GW /cm2. 

Approximately  20  percent  of  the  light  was  lost  by  reflection  from  the  in¬ 
cident  surface. 
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The  temporal  behavior  of  light  transmitted  through  the  sample 
when  damage  occurred  is  shown  in  the  multiple  exposure  cf  Fig.  31. 
Unlike  the  attenuation  observed  during  intrinsic  damage  in  ionic  solids,  10 
the  rate  of  attenuation  is  shown  for  the  pulses  which  damage  ZnSe,  and 
subsequent  laser  pulses  were  not  completely  blocked  by  the  damaged 
region.  These  observations  at  first  led  us  to  conclude  that  damage  was 
occurring  from  submicrometer  inclusions  randomly  dispersed  in  the 
focal  volume.  Such  an  interpretation  was  inconsistent,  however,  with 
calorimetric  measurements  of  the  low-field  absorption  coefficient 
(«  ~  2  \  10  cm  ).  Although  it  is  possible  that  finely  dispersed 
submicrometer  voids  may  locally  enhance  the  damage  intensity,  6^  we 
found  no  experimental  evidence  to  suggest  that  material  defects  played 
a  role  in  the  1.06  pm  measurements.  The  slow  attenuation  of  the  laser 
pulse  may  suggest  fundamental  differences  between  the  temporal  char¬ 
acteristics  of  avalanche  breakdown  in  ionic  materials  and  the  temporal 
behavior  of  intrinsic  damage  in  semiconductors.  In  the  absence  of  reli¬ 
able  dc  breakdown  data  on  ZnSe  and  intrinsic  laser  breakdown  data  at 
other  wavelengths,  it  is  not  possible  to  remove  the  ambiguity  of  our 
interpretation.  The  damage  was  well-defined,  and  any  statistical 
character  to  intrinsic  damage  was  marked  by  laser  fluctations. 

When  10.6  pm  laser  light  was  used  to  damage  ZnSe,  bulk  damage 
without  self-focusing  developed  at  Iq  =  340  MW/cm^  when  the  beam  was 
focused  to  2(Jq  =  65  pm  and  at  260  MW/rm^  when  the  laser  beam  was 
focused  to  130  pm.  Since  intrinsic  damage  limits  change  very  little  with 
laser  frequency  in  the  infrared  and  do  not  change  with  the  degree  of 
focusing,  it  appeared  that  inclusion  damage  was  occurring.  Micro¬ 
scopic  investigation  of  the  damage  sites  confirmed  this  interpretation.9'10 
It  should  be  noted  that  the  experimental  uncertainty  in  the  absolute 
calibration  of  IQ  was  ±25  percent  but  that  relative  values  of  damage 
intensities  were  determined  to  better  than  ±  10  percent.  The  inclusion 
density  p  was  later  determined  by  using  both  the  bulk  scattering  of  in¬ 
coherent  light  as  observed  under  a  microscope  and  surface  etching65 
initiated  at  inclusion  sites.  It  was  found  that  inclusion  having  diameters 
~  1  Pm  were  spaced  by  an  average  distance  of  about  45  pm  in  ZnSe. 
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Intrinsic  Laser  Damage  at  1.06ju  in  ZnSe.  Damage 
occurred  on  the  second  of  three  pulses.  The  trans¬ 
mitted  light  pulses  shown  here  are  line  drawings 
that  reproduce  the  pulse  envelopes.  High-frequency 
modulation  from  longitudinal  mode  beating  is  not 
shown. 
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Since  IQ  <1^  for  both  data  points,  the  results  of  the  preceding 

section  can  be  used  to  explain  the  damage  measurements  at  10.6  pm. 

Equation  (17)  can  be  solved  numerically  for  p  when  p  =  l/(45pm)5, 

X.  =  1.06  pm  and  =  130  pm.  -n  =  1.05  is  obtained  from  either 

Eq.  (17)  or  the  approximate  relationship  of  Eq.  (19).  Therefore, 

Ij  =  d/1.05)  (260  MW/cm^)  =  250  MW/cm^.  Inserting  the  derived 

value  of  I,  into  Eq.  (17)  with  wn  =  65  pm,  we  predict  that  InD  (65  pm)  = 

/  2  U 

330  MW/cm  as  compared  to  a  measured  value  of  340  MW/cm/cm  . 
The  experimental  data,  although  not  a  strong  test  of  the  analytical 
model  of  the  previous  section,  are  thus  seen  to  be  explained  quantita¬ 
tively  by  the  model. 

From  the  present  results  we  conclude  that  large  diameter 
10.6  pm  beams  (2uQ  ~  150  pm)  with  pulse  durations  of  140  ns  will 
induce  damage  at  250  MW/cm^.  When  other  pulse  durations  are 
used,  the  damage  intensity  will  change  if  there  exists  a  distribution 
of  inclusion  sizes.  As  noted  earlier,  IQ  oc  (tp)~2  as  long  as  Eq.  (1) 
can  be  satisfied.  An  accurate  estimate  of  the  dependence  of  1^^  on 
pulse  width  depends  on  a  knowledge  of  the  distribution  of  inclusion 
sizes  and  the  absorption  cross  section  of  each  inclusion.56  Such  a 
knowledge,  coupled  with  the  analysis  of  Hopper  and  Uhlman11  and 
the  results  of  Eq.  (17)  allow  us  to  predict,  at  least  in  principle,  the 
damage  behavior  of  a  solid  for  arbitrary  values  of  t  and  u_. 

p  0 

5.  Conclusions 

The  basic  approach  to  inclusion  damage  used  in  this  section  can 
be  useful  as  a  guide  for  materials  evaluation  studies  in  which  the  pulsed 
damage  characteristics  of  initially  transparent  solids  are  to  be  deter¬ 
mined.  By  measuring  the  intrinsic  damage  limit  I2  with  tightly  focused 
laser  beams,  a  useful  reference  is  established  for  further  damage 
studies  with  larger  focal  diameters.  Under  experimental  conditions 
such  that  damage  results  from  inclusion  absorption,  a  dependence  of 
IQ  on  the  degree  of  focusing  is  predicted  by  the  model  developed  here. 
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Careful  measurements  of  Ig  versus  Ug  and  a  determination  of  p  by 
either  light  >cattering  or  surface  etching  can  be  interpreted  in  terms 
of  the  model  of  Sec.  2  to  provide  a  characterization  of  the  inclusion 
content  of  the  solid  under  study. 
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VIII.  EFFECTS  OF  ATOMIC  IMPURITY  CONCENTRATION  ON 
INTRINSIC  BREAKDOWN  FIELDS 


Previous  measurements  of  intrinsic  breakdown  have  been  conducted 
on  pure  materials  and,  as  reported  in  Sec.  Ill,  on  one  alloy  system.  The 
possible  effects  of  atomic  impurities  on  intrinsic  breakdown  limits  have 
not  been  previously  investigated. 

In  this  section  we  report  the  results  of  a  short  series  of  measure¬ 
ments  in  which  three  laser  host  materials  both  pure  and  doped  were  damaged 
at  1.06  pm.  As  in  the  experiment  of  the  previous  section,  care  was  taken 
to  insure  that  measurable  self-focusing  effects  were  absent  and  that  only 
data  from  an  intrinsic  damage  mechanism  were  used. 

YAG  and  YAIO^  doped  with  Nd  at  concentrations  of  up  to  2  percent 
and  ruby  with  Cr  concentrations  of  up  to  0.5  percent  were  damaged  with 
the  Q-switched  laser  pulses.  The  damage  intensities  were  compared  to 
the  intrinsic  damage  levels  of  the  pure  hosts.  It  was  found  that  within 
experimental  uncertainty,  the  doped  materials  had  the  same  intrinsic 
optical  damage  intensities  as  the  respective  pure  hosts.  Incident  damage 
intensities  of  64,  71  and  210  GW/cnv  were  measured,  respectively,  for 
YAG,  YA103  and  ruby. 

+  3 

The  results  for  the  Nd  doping  were  expected  since  the  4f  valence 

electrons  of  the  trivalent  rare  earths  overlap  very  little  with  the  electron 

wave  functions  of  the  host  solid.  Avalanche  breakdown  is  a  bulk  phenomenon, 

+  3 

and  the  narrow  bands  of  Nd  effectively  remain  decoupled  from  the  host. 

+  3 

In  the  case  of  Cr  levels,  however,  substantial  overlap  exists 
between  the  valence  electrons  and  the  manifold  of  host  states.  The  Cr+ 3 
levels  are  correspondingly  broad,  and  it  might  be  expected  that  at  ex¬ 
tremely  high  doping  levels,  the  presence  of  the  additional  Cr+  3  states 
located  effectively  within  the  bandgap  will  affect  the  impact  ionization 
process  and  lower  the  breakdown  field.  On  the  basis  of  the  present  re- 
suits,  it  can  only  be  concluded  that  a  Cr  doping  level  of  0.  5  percent, 
a  value  about  10  times  the  normal  doping  level  of  pink  ruby  laser  mate¬ 
rial,  produces  a  negligible  perturbation  of  the  intrinsic  breakdown  process. 
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The  present  damage  results  along  with  the  results  of  the  previous 
section  serve  to  underscore  the  basic  insensitivity  of  electron  avalanche 
breakdown  to  material  parameters.  Although  we  have  considered  only 
two  restricted  types  of  impurities,  it  appears  safe  to  predict  that  atomic 
contaminants  of  types  other  than  those  considered  here  will  also  have  no 
measurable  effect  on  intrinsic  damage  levels. 


IX.  THEORY  OF  ELECTRON  AVALANCHE  BREAKDOWN 

A.  Introduction 

A  complete  theory  of  avalanche  breakdown  has  not  been  developed. 
The  models  of  breakdown  formulated  to  date,  which  are  based  on  relaxation 
time  approximations,  avoid  the  enormous  complexities  of  the  various 
physical  interactions  involved  in  the  avalanche  process.  Except  for  the 
work  of  Feynman  et  al.  ,  in  fact,  theories  of  hot  electron  effects  in  solids 
have  been  based  on  a  perturbation  theory  formalism  that  is,  in  principle, 
inappropriate  for  electron  avalanche  theories. 

It  is  now  recognized  that  a  realistic  theory  of  avalanche  breakdown 
in  solids  must  contain  both  a  modeling  of  the  interactions  between  the 
electron  and  the  electric  and  phonon  fields  and  a  statistical  framework  that 
accounts  for  the  distribution  of  electron  energies.  Pioneering  work  by 
Seitz  m  breakdown  in  insulators  and  Baraff  in  studies  of  electron  avalan¬ 
ching  in  semiconductors  have  shown  that  the  development  of  an  electron 
avalanche  depends  sensitively  upon  the  energy  dependence  of  the  electron- 
lattice  interaction.  A  statistical  framework  that  accurately  describes  the 
electron  distribution  function  is  thus  a  critical  aspect  of  a  relevant  theo¬ 
retical  model.  On  the  other  hand,  the  details  of  the  electron  dynamics, 
which  involve  the  quantization  of  the  EM  fields  or  of  electron  energy,  do 
not  appear  to  strongly  affect  the  predictions  of  theory.  For  this  reason, 
avalanche  models  that  use  a  classical  relaxation-time  approach  and  a 
Boltzmann  or  Fokker-Planck  equation  analysis  are  surprisingly  successful 
in  predicting  many  of  the  experimentally  observed  features  of  avalanche 
breakdown.  Such  models  are  not  rigorous  theories,  however. 

The  analysis  described  in  this  section  is  an  extension  of  a  classical 
Fokker-Planck  equation  model  of  avalanche  breakdown  that  was  originally 
developed  and  published  by  Holway 7 2 This  work  addresses,  in  a  semi- 
quantitative  manner,  some  of  the  ideas  concerning  avalanche  breakdown 
in  insulators  that  have  been  raised  by  Seitz,  FrShlich,  and  other  workers 
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in  the  field.  Among  the  ideas  to  be  tested  are  the  importance  of  the 
nonpolar  electron-phonon  interaction,  the  Frohlich  criterion  for  break¬ 
down,  and  the  sensitivity  of  breakdown  thresholds  to  the  electron  ioniza¬ 
tion  cross  section.  The  predictions  of  the  model  are  compared  to  the 
results  of  our  recent  experimental  studies  on  intrinsic  laser  damage. 
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B.  ELECTRON  AVALANCHE  BREAKDOWN  GENERATED  BY  LASER 
RADIATION  IN  INSULATING  CRYSTAL* 

Lowell  H.  Holway,  Jr.  and  David  W.  Fradinf 
Raytheon  Research  Division 
Waltham,  Mass.  02154 

ABSTRACT 

Recent  experimental  investigations  into  optical  damage  in  the  alkali 

halides  have  suggested  that  the  damage  is  due  to  avalanche  breakdown  in  its 

dc  limit.  Numerical  computations  have  been  carried  out  which  indicate  that 

the  ionization  rate  can  be  accurately  computed  by  a  diffusion  or  Fokker- 

Planck  approximation  to  the  Boltzmann  equation  for  ionization  rates  on  the 
12-1 

order  of  10  sec  .  However,  for  hot  electrons  it  is  necessary  to  add  a 
contribution  to  the  collision  frequency  which  is  due  to  the  deformation  po¬ 
tential  in  order  to  explain  the  dc  character  of  breakdown  by  1.06  pm  radia¬ 
tion. 
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1 .  Introduction 


In  a  recent  series  of  experiments  in  which  optical  damage  was  studied 
in  the  alkali  halides,  ^  v/ap  determined  that  intrin¬ 

sic  optical  damage  in  transparent  solids  occurs  as  the  result  of  electron 
avalanche  breakdown.  The  experiments,  which  used  laser  wavelengths  of 


10.6,  1.06  and  0.69  (J.m,  were  designed  to  keep  the  laser  power  below  the 
critical  power  for  self-focusing^’  ^  since,  without  this  precaution,  the 
observed  breakdown  fields  would  have  been  measures  of  this  critical  power 
instead  of  the  desired  thresholds  for  electron  cascade  breakdown.  Micro¬ 
scopic  inspection  of  the  damaged  volume  was  used  to  distinguish  intrinsic 
breakdown  from  breakdown  induced  by  inclusions  in  the  crystals.  The  ex¬ 
perimental  data  shows  that,  within  experimental  error,  laser-induced  break¬ 
down  and  dc  dielectric  breakdown  in  insulators  appear  to  be  virtually  identi¬ 
cal  processes,  a  fact  not  previously  anticipated  in  the  literature. 

Although  dc  dielectric  breakdown  in  insulators  has  been  studied  for 
nearly  half  a  century,  uncertainties  caused  by  surface  conditions,  electrode 
effects  and  space  charges  often  confused  the  experimental  results.  The  in¬ 
formation  attained  by  these  studies  was  limited  not  only  by  experimental 
difficulties  but  also  by  the  basic  inability  of  experiments  to  provide  informa¬ 
tion  other  than  dc  threshold  fields.  A  true  theoretical  understanding  of  di¬ 
electric  breakdown  from  electron  avalanching  was  inhibited  by  the  relative 
paucity  of  reliable  experimental  data  at  dc  and  the  enormous  complexity  of 
the  process  itself.  The  basis  of  this  complexity  is  the  fact  that  the  relevant 
interactions  in  avalanche  breakdown  cannot  be  treated  separately.  When 
considering  avalanche  breakdown,  it  is  not  strictly  correct,  for  example, 
to  treat  electron  excitations  as  being  distinct  from  those  of  the  lattice.  The 

electron-phonon  interaction  during  breakdown  in  NaCl  is  so  strong  that  the 

15-1 

effective  collision  rate  appears  to  be  more  than  10  sec  .  In  addition, 
distortion  from  the  ultrahigh  electric  fields  characteristic  of  breakdown 
thresholds  cannot  properly  be  treated  as  a  perturbation  on  the  energy  spec¬ 
trum  of  the  electron-lattice  system.  A  number  of  theories  have  been  formu¬ 
lated  for  dielectric  breakdown  which  largely  ignore  these  difficulties.  These 
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theories  are  based  on  an  electron  avalanche  mechanism  and  normally  use  a 
relaxation -time  approximation  to  model  the  electron-phonon  interaction. 
Such  theories  are  basically  similar  to  classical  treatments  of  breakdown  in 
gases  and  semi-conductors.  The  details  of  electron  dynamics  in  insula¬ 
tors  subject  to  ultrahigh  electrical  fields  have  been  calculated  quantum 
mechanically  by  Feynman,  et  al.,  71  but  their  treatment  in  its  present  form 
is  not  directly  applicable  to  dielectric  breakdown. 

Interest  in  dielectric  breakdown  has  recently  increased  with  the  dis¬ 
covery  that  intrinsic  optical  damage  from  laser  irradiation  develops  from  a 
mechanism  that  is  fundamentally  the  same  as  dc  breakdown.  By  utilizing 
the  great  versatility  of  the  laser,  it  has  been  possible  to  explore  the  depend¬ 
ence  of  breakdown  fields  on  frequency,  laser  pulse  duration,  and  materiai 
characteristics.  Such  studies  have  provided  new  information  on  electron 
avalanche  breakdown  that  serves  both  as  a  basis  and  a  test  for  breakdown 
theories. 

In  the  present  work  the  results  of  dielectric  breakdown  measurements 
at  dc  and  optical  frequencies  will  be  summarized  and  the  results  interpreted 
in  terms  of  electron  avalanche  breakdown.  A  Fokker-Planck  model,  which 
was  derived  without  reference  to  the  Boltzmann  equation  although  it  proves 
to  be  identical  to  the  Legendre  expansion  of  the  Boltzmann  equation70,  will 
be  used  to  provide  a  basis  for  this  interpretation.  Such  a  formalism  was 
applied  in  order  to  answer  the  following  questions: 

1.  Can  a  classical  model  of  electron  avalanche  breakdown  in  insula¬ 
tors  explain  the  basic  features  of  the  experimental  results  obtained  at  dc 
and  optical  frequencies? 

2.  How  important  are  the  details  of  the  electron-phonon  interaction 
in  determining  breakdown  fields?  In  particular,  does  the  introduction  of 
nonpolar  interactions  change  the  predictions  of  breakdown  theories  for  in¬ 
sulators  as  Seitz  suggested  many  years  ago? 
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3.  Is  breakdown  an  average  electron  phenomenon  and  /or  does  the 
avalanche  process  depend  upon  the  diffusion  of  hot  electrons  into  the  high 
energy  tail  of  the  distribution  function? 

Such  questions  cannot  all  be  answered  definitively  by  a  classical  cal¬ 
culation  based  on  perturbation  modeling  since  hu  becomes  comparable  to 
1  eV  for  the  higher  laser  frequencies  and  the  collision  frequency  will  ap¬ 
proach  10  sec'  in  the  materials  of  interest.  It  will  be  seen,  however, 
that  avalanche  breakdown  appears  to  have  the  same  fundamental  characteris¬ 
tics  whether  it  develops  in  solids,  liquids  or  gases.  For  this  reason,  sim¬ 
plified  models  of  the  process  can  be  successfully  applied.  Because  electric 
breakdown  is  similar  in  different  materials,  the  formalism  applied  here  to 
electric  breakdown  in  insulators  is  also  applicable  to  avalanche  breakdown 
in  semiconductors  and,  with  some  modification,  to  avalanche  breakdown  in 
gases,  which  is  important,  not  only  for  breakdown  caused  by  laser  beam 
propagating  through  the  atmosphere,  but  also  for  microwave  breakdown 
which  can  occur  inside  high  power  radar  waveguides. 

This  paper  begins  with  a  brief  review  of  some  of  the  important  data 
obtained  in  the  laser  experiments.  The  major  implications  of  these  data 
for  a  theoretical  understanding  of  avalanche  breakdown  are  noted. 

In  Sec.  3  and  4,  the  Fokker-Planck  model  for  avalanche  breakdown  is  ap 
plied  to  sodium  chloride-  Although  the  mathematical  formalism  for  this 
model  has  been  previously  developed  in  the  literature,  the  present  calcula¬ 
tion  differs  in  a  fundamental  manner  from  the  previous  treatment.  Where¬ 
as  the  previous  treatment  considered  only  the  polar  or  Frohlich  model  for 
the  electron-phonon  interaction,  the  present  calculation  adds  a  nonpolar, 
deformation  potential  effect  to  the  interaction  in  an  attempt  to  provide  closer 
agreement  between  theory  and  experiment.  The  calculation  is  based  on  the 
evaluation  of  the  electron  distribution  function  f(u)  and  its  growth  in  high 
electric  fields,  as  described  in  an  ionization  rate  r\  (E).  A  simplified  band 
structure  is  assumed  in  which  conduct  ion -band  electrons  can  produce  addi¬ 
tional  electrons  when  their  kinetic  energy,  relative  to  the  band  edge,  reaches 
values  that  exceed  the  band  gap. 
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Numerical  results  for  sodium  chloride  crystals  are  presented  in  Sec.  4-5. 
These  results  were  obtained  first  using  just  the  Frohlich  model  for  the  polar 
contribution  to  the  electron-phonon  collision  frequency  v  and  then  using  both 
the  polar  and  the  nonpolar  contributions  to  v  .  The  latter  contribution  is 
modeled  after  a  calculation  by  Seitz.  In  both  cases  comparison  is  made  be¬ 
tween  calculated  values  of  ti(E)  and  values  of  r|  (e)  inferred  from  optical 
damage  experiments.  It  is  found  that  the  Frohlich  model  alone  is  incapable 
of  explaining  the  important  experimental  observation  that  the  breakdown 
fields  in  NaCl  and  hence  ii  (E)  are  the  same  within  experimental  error  at 
10.6^-  and  1.06p.  On  the  other  hand,  by  adding  the  nonpolar  contribution  to 
v  and  assuming  peak  values  of  v  in  excess  of  10 15  sec-1,  reasonably  good 
agreement  is  obtained  between  the  predictions  of  theory  and  the  results  of 
experiments. 

The  present  work  demonstrates  that  knowledge  of  averag*  electron  dy¬ 
namics  must  be  supplemented  by  the  determination  of  a  distribution  function 
in  order  to  obtain  a  quantitative  description  of  electron  avalanche  breakdown. 
Avalanche  breakdown  is  a  statistical  process  in  which  the  properties  of  ele- 
trons  over  a  broad  range  of  energies  play  a  role  in  the  avalanche  develop¬ 
ment.  This  distribution,  which  is  shown  to  be  non -Maxwellian,  leads  to  a 
diffusion  of  electrons  to  high  energies.  The  diffusion  contribution  to  electron 
heating  and  subsequent  ionization  is  significant  and  can  only  be  calculated 
from  a  knowledge  of  the  electron  distribution  function.  As  a  result,  a  sta¬ 
tistical  description  of  the  avalanche  is  required  in  which  f(u)  is  calculated 
from  a  kinetic  equation.  The  present  work  suggests,  in  fact,  that  the  basic 
kinetic  modeling  may  have  as  much  effect  on  the  predictions  of  avalanche 
theory  than  as  the  precise  theoretical  modeling  of  the  electron  dynamics. 

Considerable  attention  is  given  both  to  the  physical  meaning  of  the  co¬ 
efficients  in  the  Fokker-Planck  equation  and  to  the  range  of  validity  of  its 
use  in  Sec.  4-6  and  the  Appendix.  The  relationship  between  the  present  calcu¬ 
lation  and  the  works  of  Wolf,  Shockley  and  others  will  be  discussed.  Also 
in  the  Appendix,  an  addition  to  the  kinetic  equation  that  was  ignored  in  the 
present  calculations  although  it  becomes  important  for  weak  fields,  is  des¬ 
cribed,  and  its  effect  on  the  computed  results  is  considered. 
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2. 


Experimental  studies  of  optical  breakdown 


The  experiments  invoking  laser-induced  breakdown  in  transparent 
solids  have  been  described  in  several  publications.  9’  10,  ]  1  >  18>  2°.  25,  30,  66,  67 

The  observed  damage  is  interpreted  as  an  electron  cascade  process  in  which 
the  electron  density  increases  exponentially  until  the  energy  absorbed  from 
the  beam  becomes  large  enough  to  cause  a  thermally-induced  fracture  or 
phase  change  which  permanently  damages  the  material. 


The  experiments  show  little  dependence  upon  the  laser  frequency  and 

a  characteristic  time  dependence,  which  indicates  that  the  damage  process 

is  electron  avalanche  breakdown  rather  than  multiphoton  absorption.  In  an 

avalanche  breakdown  the  number  of  conduction-band  electrons  increases  with 
time  as 


n(t)  =  Nq  exp 


J 


*1  (EXt'l  )  dt’ 


(1) 


where  E  is  the  electric  field  and  n  is  the  ionization  rate.  Figure  32  shows 
ruby  laser  light  transmitted  through  an  NaCI  sample  is  suddenly  cut  off  when 
damage  develops.  The  damage  can  occur  before,  at,  or  after  the  peak  of 
the  laser  pulse.  The  variation  between  apparently  identical  laser  pulses 
shows  the  statistical  fluctuations  in  the  breakdown  time.  The  breakdown 
after  the  peak  intensity  is  consistent  with  a  cascade  process  such  as  that 
described  by  Eq.  (1).  Since  breakdown  occurs  with  the  electric  field  within 
90  percent  of  its  peak  value,  q(E)  must  be  a  rapidly  increasing  function  of 
electric  field.  The  suddenness  of  the  attenuation  also  suggests  that  n  must 
be  large  on  a  nanosecond  time  scale. 


The  most  striking  indication  that  optical  damage  is  a  manifestation  of 
avalanche  breakdown  comes  from  a  comparative  study  of  the  electric  fields 
required  to  produce  damage  by  dc  breakdown  in  different  alkali  halides,  and 
by  comparing  different  frequency  lasers  with  10.6,  1.06  and  0.69  pm  wave¬ 
lengths.  The  breakdown  fields  for  different  materials  are  shown  in  Fig.  33.  10 
These  fields  have  been  normalized  to  the  NaCI  breakdown  fields,  which  are 
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’it?*  32  Ruby  Laser  Pulses  Transmitted  Through  a  NaCl  Sample. 

(a)  Attenuation  occurs  when  peak  laser  field  is  reached. 

(b)  Attenuation  before  peak  field  is  readied  Edamaee  = 
0-896  Epeak-  Attenuation  after  peak  laser  fielcTwas 
reached,  Edamage  =  0.954  Epeak-  (d)  Three  successive 
pulses  witli  same  peak  power,  no  damage. 
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RMS  Electric  Fields  Necessary  to  Damage  Nine 
Alkali  Halides,  Normalized  to  the  Damaging 
Field  for  NaCl. 


given  at  the  appropriate  frequency  in  Table  IV,  in  order  to  make  it  easier  to 
see  the  variation  in  trends  as  the  laser  frequency  is  increased.  The  break¬ 
down  fields  shown  in  Fig.  33  are  the  same  as  dc  breakdown  within  experi¬ 
mental  uncertainty  for  laser  wavelengths  of  1.06  pm  or  greater.  However, 
when  the  wavelength  is  decreased  to  0-69  pm,  there  is  an  increase  in  the 
breakdown  field  (for  all  of  the  halides  except  NaCl)  as  would  be  expected  from 
theories  of  frequency-dependent  avalanche  breakdown.  70  Although  the  simi¬ 
larity  to  dc  breakdown  makes  the  identification  of  the  damage  mechanism  as 
an  electron  cascade  seem  quite  certain,  the  lack  of  frequency  dependence 
for  the  1.06  pm  radiation  seems  to  require  the  hot  electrons  to  have  a  colli- 

sion  frequency  on  the  order  of  2  X  1015  sec’1  for  all  energies  up  to  the  ioni- 
zation  energy. 

-o 

Another  sequence  of  measurements,  which  indicated  that  the  ionization 
rate  h  as  a  function  of  field  is  consistent  with  dc  measurements,  69’  70  is 
compared  with  calculations  in  Sec.  5  . 

3.  The  physical  model 


The  existence  of  a  frequency  dependence  to  avalanche  breakdown  in¬ 
tensities  can  be  understood  by  a  simplified  argument.  In  the  absence  of  elec¬ 
tron  losses  from  either  defect  trapping  or  electron  diffusion,  the  essential 
features  of  the  avalanche  are  determined  by  the  rate  at  which  energy  is  fed 
into  the  electron  population.  This  rate  is  approximately  given  by  the  well- 
known  formula  for  Joule  heating 

~  =  (ne2r/m)  E2/(l  +  u2  T2)  (2) 

where  n  is  the  time -dependent  electron  density,  u  the  angular  frequency  of 
the  ac  electric  field,  and  r  the  characteristic  relaxation-time  determined 
principally  from  electron  collisions  with  lattice  vibrations.  In  general,  r 
must  be  considered  to  be  a  function  of  energy  and,  in  some  sense,  Eq.  (2) 
must  be  "averaged"  over  the  electron  distribution  function.  For  the  pur¬ 
poses  of  the  present  qualitative  discussion,  t  will  be  taken  as  a  constant, 
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TABLE  IV 


ABSOLUTE  BREAKDOWN  STRENGTH  OF  NaCl67 


Epeak  (dc) 

1.5  X  106  V/cm 

Erms  (10-6  |im) 

(1.  95  ±  0.  20)  X  106  V/cm 

E  ( 1. 06  |J.m) 

rms  ^ 

(2.  3  ±  0.  46)  X  106  V /cm 

E  (0.69  (Jim) 

rms  ^ 

(2.2  ±  0.44)  X  106  V/cm 

Cq.  (2)  shows  that  the  energy  input  to  the  electrons  scales  with  fre¬ 
quency  and  field  as  1ST /(l  .-0.  Since  the  energy  input  determines  the 

ectron  diatrtbut.on  function,  the  threshold  for  damage  will  scale  in  the 
same  way.  This  scaling  underlies  the  use  of  root-n.ea, -square  fields  in 
comparing  dc  and  laser  dielectric  breakdown.  1,  indicates  that  the  avalanche 

”7“  “7;  *°  *  *"  -  <-  «»“  'Of  laser  frequencies 

>,  <  1  tor  all  electron  energies  involved  in  the  electron  cascade.  Since 

ic  vlTuT"?  "  ‘hal  laS°r  b"l‘akd°Wn  »*  *•»«  -m  -  the  alkali  hal.de 
,  virtually  .nd.st.ng.nshable  from  dc  breakdown,  but  that  a  frequency  depend¬ 
ence  begin  to  appear  at  0.  30  micron,  it  appears  that  u  T  .  ,  at  the  ruby  la^ 
equency.  .  or  the  electron  energies  characteristic  of  avalanche  breakdown 
.erelo.e,  the  elective  value  of  r  is  about  0.  5  X  m’5  sec  or  less  for  the 

;aka‘l  '“,WeS-  ™S  effeCti-  ^“on  time  is  n.uch  shorter  than  relat¬ 
ion  tunes  interred  from  low-field  measurements  of  various  transport  prop- 
erties  in  the  same  materials.  1  P 

A  rigorous  model  of  avalanche  breakdown  would  replace  Eq.  (2)  with 
a  formalism  that  treats  the  strong  electron-lattice  interaction  quantum 

ZT 1  7  “  COUPled  eleCtr°n'lalliCe  station  easing  instead  of 

leads"  to  I  <2?  7  SlatSS  POrtUrbed  ^  ,be  latUCe-  ,he  approximation  that 
ads  to  Lq.  (2).  Except  for  the  work  of  Feynman,  et  al,  no  such  rigorous 

model  has  been  formulated  for  insulators  subject  to  high  electric  fields. 

n  o, innately,  l  ey,,, nan's  rigorous  model  of  electron  dynamics  is  not  dir¬ 
ectly  relevant  to  the  avalanche  breakdown  problem  for  two  reasons.  First 
a  statistical  framework  that  takes  account  of  the  wide  distribution  of  elec-  ’ 
ron  energies  was  not  included  in  the  calculation.  As  should  become  clear 
e  ow,  t  is  statistical  framework  is  fundamental  to  a  proper  description  of 
va  anche  breakdown.  Second,  the  particular  electron-phonon  interaction 
important  for  Feynman's  work:  namely,  the  Frohlich  interaction,  describes 
on  y  part  of  the  relevant  interaction  for  the  hot  electrons  involved  in  the 
avalanche  development. 


The  present  calculation  bypasses  the  enormous  complexities  involved 
m  the  electron  dynamics  by  using  a.,  energy-dependent  relaxation  time  model. 


A  significant  departure  from  previous  treatments  of  avalanche  breakdown 
in  insulators  is  the  introduction  of  a  nonpolar  electron-phonon  interaction. 

We  find  that  calculations  based  on  the  Frohlich  interaction  alone  are  not 
capable  of  explaining  the  important  experimental  finding  that  the  character¬ 
istics  of  laser-induced  breakdown  in  the  alkali  halides  do  not  change  meas¬ 
urably  with  laser  wavelength  for  wavelengths  as  short  as  1.06  micron. 
Following  a  suggestion  made  by  Seitz  almost  thirty  years  ago,  but  never 
tested  with  a  statistical  calculation,  we  show  that  the  addition  of  a  nonpolar 
interaction  does  explain  this  experimental  fact.  In  addition,  the  use  of 
simple  perturbation  formulas  for  the  polar  and  nonpolar  contributions  to 
the  electron-phonon  collision  frequency  v  (where  v  =  1/t)  led  to  reason¬ 
able  quantitative  agreement  with  the  measured  breakdown  fields. 

The  nonpolar  contribution  to  the  collision  frequency  is  taken  to  vary 

with  electron  energy  in  the  manner  predicted  by  Seitz  from  a  deformation 

15  - 1 

potential  calculation.  Its  peak  value  is  adjusted  to  3  X  10  sec  .  This 
contribution  is  important  to  avalanche  breakdown  in  insulators,  because 
it  is  large  at  high  electron  energies  where  the  Frohlich,  longitudinal-optical- 
phonon  interaction  is  much  smaller  than  its  peak  which  occurs  near  the  LO 
phonon  energy.  In  low-field  transport  processes  in  insulators,  very  few 
electrons  reach  such  high  energies  so  that  the  Frohlich  interaction  suffi¬ 
ciently  describes  the  electron  dynamics  for  these  processes.  During  ava¬ 
lanche  breakdown,  on  the  other  hand,  electrons  must  become  extremely 
energetic  in  order  to  gain  sufficient  energy  to  cause  impact  ionization  and 
concommitant  avalanche  growth.  The  dynamics  of  hot  electrons  thus  be¬ 
come  fundamental  to  a  description  of  avalanche  breakdown.  These  dynamics 
should  be  strongly  affected  by  the  nonpolar  interaction. 

In  order  to  make  a  clear  comparison  between  the  polar  and  nonpolar 
interactions,  we  first  carry  out  a  calculation  using  only  the  polar  interac¬ 
tion  on  defining  the  collision  frequency.  Next  we  carry  out  a  calculation 
in  which  the  deformation  potential  contribution  to  the  collision  frequency 
is  added  to  the  polar  contribution.  Only  the  gross  features  of  those  dy¬ 
namics  appear  to  be  important  and  these  gross  features  should  be  reason¬ 
ably  approximated  by  adding  the  contributions  of  Seitz'  nonpolar  interaction 
to  the  Frohlich,  polar  interaction. 
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e„ej:;::r  "i0n  T  a  SimPle  tW°-band  *>-  *•  eiectron 

SJ  pectrum.  Impact  ionization  occurs  when  a  *• 

has  gained  sufficient  energy  above  the  band  minimum  to  e'  T 
electron  across  the  band  gap.  All  en PPI,  mmimUm  to  excite  a  valence 
band  minimum.  ®aP'  AI1  ener«ies  are  measured  relative  to  the 

The  kinetic  equation 

Let  f(u’ t)du  be  the  probable  number  of  electrons  , 

tween  energy  u  and  u  4-  Hu  -+v,  ^  r  Und:  v°tume)  be- 

r°»  -  ihe  - 

tion  can  be  written  as  '  the  kinetlc  eQua- 


41  =  d  <  Au  > 

5f  7Tu~  "  t  + 

At 


1 _  /  <  (Au)2>  . 

a. .2  - - - f 


rtron:  “  *■»  - 

y  exciting  phonons  (or  mofecules)  we  have  shown78 


<  A  u> 


=  A  ( u  )  y  (u  )  -  (  hu  \ 
at' 


<(Au r>  _  4 
aT  ~  3 


=  3  u  A 


where 


2_,2 

A  =  —  E 
m 


(v  (u)  r  +  w 


Y  =  1  +2  ul2  -  v  2  u  .  9v 

3  ~  T“  (Tr)  -jt- 

w  +  v  *  v  3  u 
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and  ^  ,  the  rate  at  which  an  electron  of  energy  u  loses  energy  to  the 

dt  59  Y4  75 

lattice,  has  been  calculated  theoretically  for  alkali  halides.  ’  ’  In 

principle,  the  energy  loss  should  contribute  to  the  dispersion  also,  but  we 
have  argued  that  this  contribution  is  likely  to  be  insignificant  compared  to  the 
field  term  in  Eq.  (4b),  if  the  field  is  strong  enough,  as  is  the  case  for  elas- 
tic  energy  loss  in  electron-molecule  collisions.  ’  (When  the  field  vani¬ 
shes,  the  dispersion  must  balance  the  average  energy  loss  in  determining  a 
thermal  equilibrium  distribution,  but  the  thermal  spread  is  tiny  compared 
to  the  spread  when  breakdown  is  caused  by  a  strong  field.) 


It  can  be  shown  that  the  field  dependent  terms  in  Eq.  (3)  are  identical 
with  the  terms  which  are  obtained  by  expanding  the  Boltzmann  equation  or 
Legendre  polynomials70,  77,78,  when  f  is  considered  a  function  of  time. 

One  of  the  advantages  of  the  Fokker- Planck  derivation,  besides  some  reduc¬ 
tion  in  formal  manipulation,  is  that  the  conditions  for  its  applicability  are 
well  known.  For  example,  it  is  necessary  that  the  energy  gained  in  one  col¬ 
lision  time,  v  "1,  be  much  less  than  the  total  energy,  which  amounts  to  the 
condition 


e2E2 

m(v  +  u  ) 


«  I 


(7) 


where  I  is  the  ionization  energy.  In  other  words,  the  diffusion  formula¬ 
tion  is  limited  to  electric  fields  which  are  weak  in  the  sense  of  Eq.  (7). 

It  is  not  difficult  to  show  that  if  Eq.  (7)  is  violated,  the  electron  velocity 
will  take  on  a  considerable  component  antiparallel  to  the  field  in  the  time 
between  collisions  so  the  velocity  distribution  function  will  no  longer  be 
close  to  spherically  symmetric.  Thus  the  Legendre  expansion  method 
must  similarly  be  expected  to  fail  for  strong  fields,  which  is  consistent 

with  our  assertion  that  the  Legendre  expansion  and  the  Fokker-Planck  ap- 

2  2  30 

proximation  are  identical.  If  we  take  1=7.7  eV  and  v  +  u  =10  ,  the 

equality  in  Eq.  (7)  will  only  occur  for  fields  as  large  as  66  MV/cm,  which 
indicates  that  the  Fokker-Planck  equation  is  valid  for  the  measured  break¬ 
down  fields. 
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The  amount  of  energy  per  unit  volume  absorbed  by  the  electrons  can 
conveniently  be  obtained  by  multiplying  Eq.  (3)  by  u  and  integrating  over 
u.  The  second  term  on  the  right  hand  side  vanishes  upon  integrating  by 
parts,  leaving 

w  \  ufdu  =  j  fdu 

o  o 

oo 

=  ^  A(u)  Y (d)  fdu 
o 

where  we  have  dropped  the  energy  loss  term  in  Eq.  (4a).  Equation  (8)  denotes 
the  total  energy  gained  per  unit  volume  from  the  field,  and  can  be 
written  as 


9  - 

3T  nu 


2  _2 
e  E 

m 


n  v 


eff 


2  2 
veff+  U 


(9) 


where  is  not  the  average  collision  frequency  (except  in  the  trivial 

case  v=  constant).  For  example,  if  v  =  bu  and  if  f  were  a  Maxwellian 

5 

at  temperature  T,  we  would  obtain  bkT  for  u  »  and 

veff  =  "2  bkT  f°r  u  <<:  veff"  This  last  result  is  known  in  ionospheric  propa¬ 
gation,  where  it  is  obtained  after  considerable  algebraic  manipulation  of 

79 

the  integrals  arising  from  the  Boltzmann  equation.  In  general,  if 
v  =  bum  and  m  >  -3/2,  vnFe  =  (1  +  4  m)  v  ,  for  u  »  v,  which  can  be  shown 

0x1  O  3,  V  0 

from  Eqs.(5),  (6)  and  (8)  and  holds  for  nonequilibrium  distribution  functions. 


In  order  to  carry  out  numerical  calculations,  we  must  assume 
some  model  for  v  since  no  entirely  satisfactory  model  for  hot  electrons 
exists.  The  first  model  considered  will  be  that  of  Frohlich74  which  is 
derived  by  assuming  that  positive  and  negative  charges  situated  at  the  lattice 
sites  are  perturbed  by  the  electrons  and  has  the  form 


V,  -1  -1/2 

v  =  b  u  u  ' 

u  <  u 

o 

i  -3/  2 
v  =  b  u  ' 

u  >  u 
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for  electrons  with  energies  somewhat  greater  than  thermal  where 


(11) 


and  uq  is  a  crossover  energy,  approximately  equal  to  the  electron  energy 
corresponding  to  the  boundary  of  a  Brillouin  zone,  which  we  take  as 

UQ  =  22/3  it 2  h2/8ma2  (12) 

Our  numerical  examples  are  for  NaCl  crystals  where75  the  nearest  inter¬ 
ionic  distance  a  =  2.81  A,  the  effective  mass  is  taken  as  the  electron  mass 
m  =  9.109  X  10"28,  ut  =  3.09  X  1013,  T  =  300°K,  and  M  =  23.2  X  10-24gm 
Equations  (11)  and  (12)  are  Frohlich's  original  formulas74  which  are  larger 
by  a  factor  of  2  than  values  given  by  O'Dwyer,  a  numerical  factor  does 
not  affect  the  results  greatly.  This  collision  frequency  is  plotted  as  curve 
(1)  in  Fig.  34,  where  uq  =  1.89  eV. 

For  nonpolar  crystals,  the  deformable  atom  hypothesis  can  be  used 
to  obtain  a  perturbation  potential  and  Seitz59  indicates  this  should  be  im¬ 
portant,  even  in  polar  crystals,  for  hot  electrons. 

This  contribution  to  the  collision  frequency  has  the  form 


v  =  d  u 
s 

u  <u0  . 

V  .  du  3/2  u -1'2 

s  o 

u>uo  , 

(13) 

for  energies  much  greater  than  thermal.  For  our  purposes,  we  consider 
d  a  phenomenological  constant  and  add  the  contribution  of  Eq.  (13)  to  the 
polar  contribution  given  by  Eq.  (10),  so  that 

v  =  buQ  '  u  ^2+du  u<u 

v  =  b  u"3/2  +  d  uQ3/2  u  -1/2  u  >  u0  •  (14) 
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14 

Curves  2  and  3  are  the  result  of  choosing  duQ  =  6.  5  X10  (so 

that  (10)  and  (14)  make  equal  contributions  at  u  )  and  choosing 
15  ° 

du  =2X10  .  Curve  4  is  the  form  for  a  constant  mean  free  path  arbi- 

°  15 

trarily  normalized  to  10  at  u  =  1  eV.  Curves  1  and  4  bracket  the  extremes 

”3/2  ^  r) 

used  in  the  literature  since  the  u  '  form  (Curve  1)  was  used  by  Keldysh 

77 

while  the  constant  mean  free  path  assumption  (Curve  4)  was  made  by  Wolff 
and  Baraff.^ 


The  net  number  of  electrons  per  unit  volume  whose  energy  increases 
from  a  value  less  than  u  to  a  value  greater  than  u  per  unit  time  can  be 
defined  as  J(u),  the  current  in  energy  space.  Then  Eq.  (3)  is  equivalent  to 


Of  _  aj 

cRT  3u~ 

where 

T  _  <  £u  >  ,  1  9 

J_  f  "  2"  ciu 


<  (Au)2  > 
£t 


(15) 


=  Gf  -  d|1 


(16) 


where  G  = 
and  D(u)  = 


1 


I 


A  -  (^)l  Play?  a  role  like  an  average  rate  of  energy  change 
uA  plays  the  role  of  a  diffusion  coefficient  in  energy  space. 


Equation  (16)  is  a  convenient  form  for  numerical  purposes  since 
>  the  derivative  of  the  collision  frequency  has  cancelled.  It  is  important 
to  note  that  G  is  not  the  rate  of  change  of  energy  for  the  average  electron 
which  is  given  by  Eq.  (3)  in  general  (and  equals  A  -  (^j-)L  if  v  is  constant). 

If  diffusion  could  be  neglected,  breakdown  would  only  occur  if  G  >  0. 
Thus  we  find  it  more  meaningful  to  plot  A  and  contrast  it  with  the  energy 
loss,  instead  of  contrasting  A  with  the  energy  loss.  This  has 

been  done  in  Fig.35  a  where  Curves  1  and  2  plot  ^  A  based  on  the  collision 
frequency  of  Curve  1  of  Fig.  34.  The  energy  loss  rate  is  given  by  Curve  3 
so  that  G  >  0  if  Curve  1  exceeds  (3).  Curves  4  and  5  plot  D(u)  =  jj-  1—^2 — — 
for  different  fields.  The  diffusion  constant,  which  has  units  eV^/psec, 
cannot  strictly  be  compared  with  G  but  its  relative  size  plays  an  important 
role  as  can  be  shown  by  the  following  qualitative  argument.  If  steady-state 
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Collision  Frequency.  Curve  (1)  is  Frohlich's  model 
from  Eq.  (10):  curve  (2)  includes  nonpolar  contributions 
of  Eq.  (14)  with  d  u0  =  0.65  X  10l5;  curve  (3)  has  d  u0  = 

2  X  1015;  curve  (4)  is  constant  mean  free  path  =  1  X  1015 
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Fig.  35a  Coefficients  in  the  Fokker-Planck  Equation  for 
Collision  Frequency  (1)  with  \  =  1.06  fim. 
Curves  (1)  and  (2)  give  A/3  and  (4)  and  (5)  give 
D(u).  Curve  3  is  the  energy  loss  (d  u/dt)T 
Curves  (1)  and  (4)  are  for  E  =  10  M\  /cm  while 
(2)  and  (5)  are  for  15  MV/cm. 
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ionization  is  going  on,  ,T  must  be  greater  than  zero  for  u  between  a 

value  Uj  near  zero  and  a  value  near  the  ionization  energy  I.  Since  I)  is 

0  f 

positive,  if  G  ,s  negative,  then  must  also  be  negative  with  a  magnitude 


9f  >  G  f 
15TT  I  D  f 


which  implies 


f  <  exp 


.  G  ! 
1  n  i 

l  u  i 


d  u' 


(17) 


u , 


Thus,  if  G  is  negative  and  -  (G/D)  I  is  much  larger  than  one,  the  number 
of  electrons  reaching  ionization  energy  will  be  exponentially  small.  Figs.  35a 
and  b  compare  the  coefficients  for  X  =  l .  06  and  X  =  l .  06  pm.  For  E  =  10 
Mv/cm  in  Fig.  35a,  we  see  from  Curves  1  and  3  that  G  is  negative  while  D 
is  less  than  -G  below  1  eV  and  is  only  slightly  larger  above  1  eV.  A  has  its 
maximum  near  0.3  eV  where  v  =  u.  On  the  other  hand,  A/3  is  much  larger 
than  (d  u/dt)j^  in  Fig.  35b  except  below  1  eV  and  has  a  maximum  near  4.  4  eV 
where  v  =  u  .  In  addition,  D  is  extremely  large  except  near  the  origin  so 
we  expect  the  ionization  rate  to  be  orders  of  magnitude  larger.  Note  that 
it  is  difficult  for  these  electrons  below  0.  5  eV  to  gain  energy,  so  there  is  an 
effect  like  a  low  energy  trap. 

The  energy  loss  term  used  in  Fig.  35  was  the  form  given  by  Frbhlich74’  75’ 
neglecting  deformation  potentials,  which  is 


7 


(3F)L  =  dl  u"l/2en 

where 

7  3  (4/fi  u  )u 
7  =  (4  u J/ 2/*  ut)  u1/2 


(18) 


u  <  u 


u  >  n 


and 


dj  =  Tie4  (2m)1/  2/  Ma3 


The  current  J  is  positive  so  that  electrons  are  carried  up  beyond  I 
the  ionization  potential,  given  as  7.7  eV  for  NaCl. 
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An  electron  having  an  energy  u,  greater  than  I  can  use  that  energy  to 
move  a  second  electron  from  the  valence  band  to  the  conduction  band  by 
impact  ionization,  leaving  an  energy  u-I  to  be  shared  by  two  electrons  in 
the  conduction  band.  If  the  energy  is  ahared  equally,  Eq.  (15)  becomes 

0  £  0  j 

=  ■  ”  w(u)  f(u)  +  4ia  (2u+  1)  f(2u  +  I)  ;  (19) 

Integrating  Eq.  (19)  over  u,  the  current  term  vanishes,  yielding 

00 

|f  =  \  V.(u)  f(u)  dy 
I 

where 

oo 

n  =  \  fdu 
o 

^  or  lack  of  information  on  v^(u)  we  assume,  not  only  that  it  is  zero  below  I, 
but  also  that  it  increases  linearly  above  I  as 


w(u)  =  a1  (u  -  I)  u  >  I 

where,  in  the  numerical  examples,  we  have  taken  a1  =  4  X1014sec_1  eV-1. 
this  value  insures  that  most  electrons  will  cause  ionization  within  1  or 
2  eV  of  I,  so  that  the  results  are  only  slightly  different  from  those  obtained 
previously  by  assuming  that  all  electrons  reaching  the  ionization  energy 
immediately  cause  ionization.  In  atmospheric  breakdown  Kroll  and  Watson70 
use  \  G(u  |  u*)  f(u* )  du*  in  place  of  the  last  term  in  Eq.  (19)  which  would 
i educe  to  our  term  if  they  took 

G(u  |  u' )  =  4  (u  1 )  6  (u*  -  2u  -  I) 
as  would  be  proper  if  the  energy  were  shared  equally. 


The  boundary  conditions  for  Eq.  (19)  are  J(0)  -  0  and  f  approaches 
zero  as  u  -  <*.  We  anticipate  that  the  solution  for  f  can  be  written  as 


a  (he 
m 
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where  <f>  and  i»m  are  eigenfunctions  and  eigenvalues  of  Eq.  (19),  i  .  e.  , 

-  J-  JU  )-v.(u )  <b  (u)  +  4v.  (2u+I )  6  (2u+I)  =  ca  <j>  (u)  ,  (29) 

9u  vm  i  mn  i  vm  m  m 

where  E  is  assumed  to  remain  constant. 


The  a  can  be  determined  from  the  initial  condition.  If  a  is  the 
m  o 

largest  eigenvalue,  at  late  times  the  distribution  function  assumes  the 
asymptotic  form 

f  =  a  (u)  exp  (a  t)  (21) 

o  o  o 

which  means  that  the  number  density  satisfies 


9n 

c)T 


T]  (E)  n  =  a  n 
o 


(22) 


We  have  written  a  finite  difference  computer  program  in  which  we  specify 
an  initial  distribution  function  and  follow  its  evolution  in  time  by  means 
of  Eq.  (19).  After  increments  At  we  calculate 


a  =  — 
At 


/  n(t) 
n(t-At) 


) 


(23) 


and  find  that  a  converges  to  a  constant  value  a  =  n  (E)  sec  1  after  times 

-1  ° 

which  are  typically  of  the  order  of  0.  2  rj  .  Since  this  is  very  fast  com¬ 
pared  to  the  time  in  which  the  electric  field  changes,  we  can  calculate 
rj  (E)  for  fixed  E  and  then  compute  the  electron  density  from 


n 


r  t 


nQ  exp 


( 


U  (E(t')  )  dt' 


as  anticipated  in  Section  2.  We  may  sometime  refer  to  Eq.  (2l)  as  a  "steady 
state"  even  though  it  is  important  to  realize  n  is  exponentially  increasing 
and  the  current  is  not  zero. 


Ill 


Numerical  results 


Numerical  results  for  NaCl  are  Siven  in  this  section.  They  differ 
rom  the  results  obtained  previously  for  sapphire72  because  v  is  about  an 
order  cl  magnitude  larger.  Thus,  for  0.69  radiation,  u  »  „  except  for  en¬ 
ergies  below  0. 1  eV.  On  the  other  hand.  Prohlichs  model  for  NaCl  has  r  -u 
near  0.3  eV  for  1.06  pm  lasers  and  v  .  u  near  4.4  eV  for  10.6  pm  lasers. 

oreover,  since  uo  is  0.8  eV  for  sapphire,  it  was  reasonable  to  use  the  col¬ 
lision  frequency  formula  defined  for  u  >  u„,  but,  for  NaCl.  un  =  1<89  eV  so 

a  we  must  also  use  the  u  <  uQ  formula.  Another  difference  from  the  elrlier 
calculation  is  that  we  now  include  an  ionization  frequency  ,,  as  discussed  in 

by  asstmin?0"’  T  giV6S  S"BhtlylOWer  -tea  than  are  obtained 

assuming  every  electron  whose  energy  reaches  I  causes  impact  ionization. 

in  whic?  I?'"  difrnCe  Pr°Sram  WaS  adVanCed  by  an  6Xplicit  CuIe'-  nlethod 
to  maint  •  t  h-i-  *  ^  tyP1Cally  as  0,1  c'’  and  At  was  kept  less  than  (Au)2/2d 

two  and T,nb  f  ty'  0,16  *“*  Pr°Wem  *  “  by  a  t actor  of 

two  and  At  by  a  factor  of  four,  which  resulted  in  1. 4  percent  change  in  the 

ca  CU  ated  value  of  „  .  The  program  was  written  with  a  second  order  llunge- 

utta  integration  as  an  option,  but  this  nearly  doubled  the  running  time  with 
only  a  slight  increase  in  accuracy. 

trie  filidr  15  SMw7  a  diS!ribUti0n  “  that  ™  calculated  for  an  elec- 

ti  '  /cm  and  a  1.06  pm  wavelength,  with  an  initial  distribu- 

assumed  to  be  Maxwellian  at  300°K,  with  an  electron  density  of  108  per 
vo  ume.  After  0.03  psecs  the  distribution  function  has  already  spread 
oonsi  era  ly  from  its  initial  shape  which  was  concentrated  around  .025  eV. 

The  evolution  of  the  distribution  function  can  be  seen  more  clearly 

r,0m  etd7ity  C°ntOUrS  diSplayGd  in  Pi^  87a  and  b.  Here  a  quantized 
representation  has  been  used  in  .  q  ° 

a  base  of  106  eV1.  -f  r  W6  Pn  an  inteRer  "°"  if  f  is  below 

,  , creeds  the  base  by  less  than  1  dB.  the  entry  is 

eft  blank  Por  1-2  dB  above  the  base  a  is  printed  and  for  2-3  dB  above 

d„  .  aSa  ?nt''y  iS  blank‘  The  alS°ri*hm  proceeds,  leaving  every  other 
dB  blank  and  printing  an  integer  "N"  when  f  exceeds  the  base  by  an  amount 
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Distribution  Function  which  Evolves  from  an  Original 
Maxwellian  at  300°K  and  a  Total  Number  Density  of 
10®  per  Unit  Volume  when  NaCl  is  Irradiated  by  an 
RMS  Field  of  15  MV/cm  with  a  1.06  M-m  Wavelength. 
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Fig.  37a  Quantized  Contours  of  a  Distribution  Function  Which 
Evolving  from  an  Initial  Maxwellian  at  300° ,  Mrhen 
NaCl  is  Irradiated  by  a  1.06  pm  Laser  with  E  =  15 
MV/cm.  The  first  column  represents  the  initial  con¬ 
dition  when  almost  all  the  electrons  are  on  the  lowest 
energy  bin  and  the  distribution  spreads  to  a  asymptotic 
form  after  about  .  14  psecs.  Collision  frequency  is 
curve  (1)  of  Fig.  34. 
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Fig.  37b  Quantized  Contours  Evolving  from  a  Uniform 
Distribution  Function,  E  =  10  MV/cm,  \  = 

1  *  06  pm.  With  a  lower  field,  the  asymptotic 
distribution  function  is  compressed  in  toward 
lower  energy. 
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ranging  from  2N-1  to  2N  dB.  Finally,  a  "9"  is  printed  when  f  exceeds  the 
base  by  more  then  17  dB.  Each  column  represents  the  distribution  function 
at  a  specific  time.  Thus,  the  16th  column  in  Fig.  37a  gives  the  same  in¬ 
formation  as  the  curve  for  .03  psecs  in  Fig.  36.  Each  horizontal  line  gives 
the  density  as  a  function  of  time  at  the  center  of  an  energy  bin  0. 1  eV  in  ex¬ 
tent.  Only  every  second  bin  has  been  printed  so  the  bottom  row  corresponds 
to  0.05  eV,  the  next  row  to  0.2  5  eV,  etc. 

Figure  37a  shows  the  evolution  of  an  initially  Maxwellian  distribution 
function  for  the  same  conditions  as  Fig.  36.  The  electrons,  irradiated  by  a 
field  E  =  15  MV/cm  which  is  switched  on  att=0,  spread  in  energy  space  until  they 
have  nearly  reached  an  asymptotic  form  after  0. 14  psecs.  Figure  37b  shows 
the  compression  of  an  initial  distribution  function  f  =  108/Iforu  <1  when 
NaCl  is  irradiated  by  a  10  MW/cm  field  with  1.06  pm  wavelength.  With 
the  weaker  field,  the  distribution  function  at  late  times  collapses  closer  to 
the  origin  than  the  asymptotic  distribution  function  in  Fig.  37a. 

The  electron  currents  in  energy  space,  normalized  to  one  electron, 
are  plotted  m  Figs.  38a  and  b  for  the  same  conditions  as  those  contained  in 
bigs.  37a  and  b.  In  Fig.  38a,  the  electrons  initially  flow  towards  higher 
energy,  and  finally  reach  the  steady  state  form  after  0. 15  psecs. 

The  current  J(I),  where  I  is  the  ionization  energy,  has  a  special 

significance  which  can  be  seen  by  integrating  Eq.  (19)  from  0  to  I,  which  pro¬ 
duces 

.  *  .1  31 

3n  (  d  f  ,  .  (  ( 

7TT  ~  p  FT  ~  “  j  v  .f  du  +  2  J  v.(u)f(u)du 

where  the  boundary  condition  J(0)  =  0  has  been  used.  The  second  term 
on  the  r.  h.s.  vanishes,  and  since 
31  oo 

j  v i  f  du  ==  j  v.fdu- 

we  obtain 


116 


E  =  10  MV /cm 


— —  Negative  flux 
—  —  Positive  flux 


i  ,  j,  IU,JIIL|„|M  iiUUM  .1 WI-MU  Ik  I 

’ 


3n 

Tt 


=  J(I) 


oo 

I 

I 


3f 

sr 


du 


(24) 


which  states  that  every  electron  passing  I  either  causes  an  ionization  or  else  is 
used  to  increase  the  number  of  electrons  in  "storage"  at  energies  greater 
than  I.  Since  the  last  term  is  relatively  small,  we  have  approximately 


J(I)  ~  1  3n 
n  n  <H" 


=  n  (t) 


(25) 


where  n  is  the  instantaneous  ionization  rate.  For  example,  from  Fig.  (38a) 
the  asymptotic  value  of  J(I)/n  at  7.7  eV  is  0.74  psec"1  which  compares  with  a 
numerical  result  for  r)  of  0.73  psec  at  0.15  psec. 


In  Fig.  38b  the  current  is  negative  at  early  times  as  the  electrons 
rapidly  lose  energy,  except  for  electrons  in  the  tail  of  the  distribution 
which  diffuse  upward  and  are  eventually  ionized.  Finally  after  1.0  psecs,  the 
distribution  reaches  its  asymptotic  form  and  J(u)  is  everywhere  positive.  Note 
that  the  asymptotic  form  of  J(u)  does  not  monotonically  decrease  as  it  would 
do  if  all  the  electrons  were  born  at  zero  energy.  Instead,  a  considerable 
amount  of  electrons  are  created  with  energies  as  high  as  1  eV  which  re¬ 
quires  J(u)  to  increase  at  low  energies,  as  shown  by  Eq.  (19),  suggesting 
that  a  considerable  fraction  of  electrons  avoid  the  low  energy  "trap"  . 
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The  calculated  ionization  rates,  based  upon  curve  1  of  Fig.  34, 

for  the  collision  rate  frequency  and  Frohlich's  model  for  the  energy  loss 

are  shown  in  Fig.  39  where  they  are  compared  with  experimental  results'1'  for 

1.06  pm  pulses.  The  experimental  ionization  rates  were  taken  equal  to 

1 8 / 1  where  t^  was  the  measured  pulse  length  for  breakdown.  (The  factor 

"18'^was  obtained  by  assuming  the  ambient  electron  density  must  be  in- 
8  17 

creased  by  10  to  cause  damage.  This  factor  would  be  increased  to  "  27" 
if  Seitz’  "  40  generation"  model  was  used.) 


< 

1 


2 

If  the  time  variable  in  Eq.  (2)  was  changed  to  t'  =  tE  ,  and 
(du/dt)^  were  neglected,  then  the  electron  field  would  disappear  from 
the  equation.  When  the  electric  field  becomes  large  enough,  this 
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Fig.  39  Ionization  Rate  from  Theory  and  Experiments  for 

1.06  and  10.6  pm  Radiation.  The  solid  curves  con¬ 
nect  points  calculated  using  curve  (1)  of  Fig.  34  and 
Frohlich' s  loss  rate.  The  error  bars  indicate  experi¬ 
mental  points  and  the  dashed  straight  lines  are  loci 
for  ionization  rates  proportional  to  E2.  Here  E  is  the  rms  field 
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neglect  will  be  valid,  thus  t)  will  become  proportional  to  E^.  This  expec¬ 
tation  is  born  out  be  Fig.  39,  where  the  dotted  straight  lines  were  drawn 
as  loci  r|  / E  =  constant  and  correspond  to  high  field  asymptotes  to  the 
theoretical  curves. 

Considering  that  there  were  no  adjustable  constants  in  calculating 
the  curves  in  Fig.  39,  the  agreement  may  be  considered  quite  good.  How¬ 
ever,  our  other  experiments67indicated  little  difference  between  1.06  |im 
and  10.6  fim  radiation, which  could  only  be  exactly  true  if  v  >  >  u  at  all 
energies.  Since  Seitz'  '  suggested  the  deformation  potential  will  be 
important  for  hot  electrons,  we  have  made  further  calculations  using 
curve  3  of  Fig.  34. 

To  be  consistent,  deformation  potential  effects  should  be  included 
in  (du/dt)^,  so  we  have  taken 

/ du ,  _  ,  du  v  ,  ,du  v 

(HT)L  "  {HF)1  +  (HT)2  (26> 

where  (du/dtlj  is  the  polar  crystal  effect  given  by  Eq.  (18)  and 


(  du, 

(HF'2 


=  (3/ 2^c  n/  2mu  vs,  u  <  uQ  , 
=  ^3/  2^c  n/  2muQ  v  g  u  >  uQ  , 


u  >  u. 


where  v  g  is  the  deformation  potential  collision  frequency  from  Eq.  (13) 

and  c  =  velocity  of  the  longitudinal  optical  phonon  =  4.58  X  105  cm/sec 
82 

for  NaCl.  Equation  (27)  follows  from  requiring  Eqs.  (10)  and  (11)  of 
5  9 

Seitz  paper  to  be  consistent  with  Eq.  (13). 


The  energy  loss  from  Eq.  (26)  is  plotted  as  curve  3  of  Fig.  40 
(which  is  similar  to  a  sketch59  given  as  Seitz'  Fig.  34b).  The  other 
curves  are  plots  of  A/3  and  D(u)  which  are  based  on  curve  3  of  Fig.  34 

c  q 

(which  is  similar  to  a  sketch  given  as  Seitz'  Fig.  33b). 


This  plot  is  for  \  =  1.06  but,  since  v  exceeds  u  over  most  of 
the  energy  range,  the  coefficient  do  not  change  drastically  for  \  =  10.6. 
The  coefficient  G  (and  the  rate  of  energy  change)  is  negative,  so  the  break¬ 
down  occurs  due  to  the  diffusion  term,  D(u). 
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The  results  in  Fig.  41  include  deformation  potential  effects  and 
were  computed  with  curve  3  of  Fig.  34  for  the  collision  frequency  and 
Eq.  (26)  for  the  energy  loss  to  the  lattice.  The  RMS  field  is  8  MW/cm  and 
the  wavelength  is  1.06  pm.  Figure  41a  shows  quantized  density  contours 
tor  an  initially  uniform  electron  density  collapsing  towards  its  asymptotic 
form.  Figure  41b  shows  the  energy  flux  J(u)  at  different  points  in  time. 

The  flux  at  .01  psecs  is  negative  below  6.2  eV  since  most  of  the  hot  elec¬ 
trons  lose  energy  initially.  At  0.68  psecs,  the  distribution  appears  to  be 
at  its  asymptotic  limit,  judging  from  the  convergence  of  successive  values 
of  r|  .  it  is  possible  that  a  steady  state  could  exist  in  which  many  electrons 
are  born  near  0.  5  eV  and  lose  energy  in  order  to  populate  the  energy  range 
below  0.5  eV,  with  f  increasing  as  <f>  g(u)ea^  .  This  may  explain  the  nega¬ 
tive  flux  below  0.4  eV  in  the  0.68  curve  in  Fig.  41b.  However,  because  we 
have  not  carried  this  calculation  to  later  times,  we  cannot  be  sure  whether 
J  becomes  positive  at  later  times. 

Figure  42  plots  the  ionization  rate  including  the  deformation  potential 
effect  in  the  collision  frequency  with  the  solid  curves  showing  results  using 
the  energy  loss  rate  given  to  Eq.  (26)  while  the  dashed  curve  uses  only  the 
Frohlich  energy  loss  from  Eq.  (18).  The  figure  shows  clearly  that  when  E 
becomes  large  enough  so  that  the  loss  rates  are  inconsequential,  the  curves 
approach  their  asymptotic  forms  proportional  to  E  .  Compared  to  Fig.  39, 
the  10.6  pm  curve  has  moved  downward  while  the  1.06  pm  curve  has  moved 
upward  so  that  the  two  curves  are  close  enough  together  so  that  the  break¬ 
down  at  1.06  pm  would  appear  to  be  dc  within  the  experimental  error.  The 
agreement  is  good  except  for  the  experimental  points  near  the  bottom  of  the 
figure,  which  have  larger  ionization  rates  than  would  be  expected  from  the 
calculations.  The  calculated  points  would  move  upward  if  (du/dt)T  were  de- 
creased,  or  if  the  loss  processes  responsible  for  (du/dt)L  cause  a  disper¬ 
sion  which  is  an  important  contribution  to  the  diffusion  coefficient  D  for  fields 
less  than  about  5  MW/ cm. 

Once  again  the  ionization  frequency  for  an  electron  with  energy  u  was 
taken  as  aj  ^u-I)  where  a^  was  4  X  10^  sec  ^  eV  This  was  large  enough 


123 


AT)  (roTTfinFriTi  oo  oTFO  0  o  o  p'p  o  o  ifnn  oPoponooTfpooPoPo'oopooonooppooor,  pooppopopooppp 
00  0_0  0  0  0  000  o  fl  POO  0  0  o  o  o  o  0  0  0  o  o  p  o  o  o  o  o  o  o  p  o  o  o  r,  p  o  o  gp  p  p  o  p  o  p  p  o  p  p  p  n  p  p  p  p  p  p  p  p  p  p  p  p  p  p 
0  0  o  o^  o  o  li  o  o  o  oo  o  oo  o  n  n  n  pop' o"o  TTo  pp-pppppporiflp  plTp  poooopoopnoooppnnPOOoopppp 
Ql  ill  oofln ooopoopppppppppppAOPnppopooAApoppooppppftpppppppppooonnpppp 

0???  TT1  Oooooopp o  P(TVoo"0"Oo7fo  p"p P pp"p 071Y  p p p"p pffp'p'ppTft  p  0  p  p  p'o'oVp'OTi  07 ffp  p  o  o 
013  ??  Ill  PPPPPPPOPOPOPPPPPPPAPOPPPPPPPPPPPPPPPPPPPPPPPPPOPPPPPP 

A-  "  33  2?"2  Ill  ^'oflO'flprtOpppp'ppppppOpppppppnpppppppppppppppppppppnpp 

6  U  U  33  ?  ?  ?  1111  PPPPPPPPPPPPPOOPPPPPPPPPPPPPPOPPPPPPPPPPPPPPPP 

b  uu  5V  ??22  ilM  0 p a  op pnpppo nnnpp p popp ppp p nn pp pp p p p p op pp pp p p 

65  u  u  3  3  3  ?  2  ?  1111  oppppppppppppppppppppppppppppppppppppppp 

655  uTu  XXXT  ?T?’?  TT'H  ‘  ppppp'ppopoppPpppppp'pppppPpp'pCPPpppppp 

5_  55 _ Uuu  _  333  ??2?  iili  p  p  OOP  pop  p  pP  ppppppp  p  p  p  p  p  pp  p  p  o  pp  p  p  p  p  p 

6  55  uuu  333  ????  1111  PPPiippppppppppppppppppppppppppppp 

6  555  uuu  3333  ????  111i  pppppoppppppppppppppppppprpnpp 

p  555  uuu  35i  ?  ?  ?  ?  ?  Iili  pppppopppppppppppppppppppppp 

64  65  4  5_ u  uu  3  3  3  3  ?  ?  ?  ?  ill!  ppooppppppppppppppppppppppp 

66  65  5  u  u  u  ,7  113  3*  ?  ?  ?  ?  11111  Pnprnpnnn?ip’p"pppnppppppnpp 

66  555  u  UUU  <333  ?  ?  ?  ?  1111  pppppppp  pppppppppppi-f' 

66  5556  U  U  u  u  3  3  3  3  ?  ?  ?  ?  ?  1  M  M  ppppppppppppppppppppp* 

66  6556  UU,I  4  3333  ????  i'll  i'PPPPP*.  pppppppfpp  ppp 

6  6  6  6566  U  4  U  ,J  <333  3  PP??  Mill  opppppppppppnrppnpp 

6  6  6  56444  U  6  ,1  4  3  3  3  3  3  ?????  11(11  PPPPPPPPPPPPPPPPPP 


66 
646 
4  6  6 

7T5 — 
666 
666 
666 
66 
66 

6T? - 

6  6 


556  6 

it  l;UU 

64446  uii'iuu 
56  456  11  JUUU 

46455  ju;| 

446556  uu 

65566  uu 

■  “"555555  u 

6466444 
455556 
5654544 
“  5565644 

^  4444444 

■■•‘41,644546 


uu  U  u 
U  U  U  4  4 
UU4  44 
UUU  U  U 

li  uuu  » 


33333  ?P???  |  1  '1  1  1  PPPPPPP  0  PPPPPpP  A 

33333  P  ?  ?  ?  ?  1111  pPPPpuPPPipPpp  .* 

33333  PPPP?  1111  uPPPPPPf,  PPP6PP 

33333  ?P?P  11111  P*PPPfPfPrpr 

33333  PPP?  Mill  P  0  0  0  0  f,  0  0  n  0  0 

J  3  3  333  ????  1  1  1  1  1  PPiPPPPPppP 

14  "3333  3  ??P?  11  Ml  OPPPPPppP 

UUU  3333}  ?  ?  2  ?  11111  ('PfPPPPP 

3  3  3  3  3 


UUOU.'J  33  333  ?????  Ill  11  Pppppp 

UUUU'14  3  3  3  33  ?????  1)111  PPP, IP 

U  U  ,1(1  -J  3  3  3  }}  ??????  11M1  PPP 

_ 5f_UUU  3  3  3  3J3  p?????  Mill 

"IfuuuuV  '  3  33  3  3  2222?  11111 


66  45446444  ""Tfuuu'iV"  "  3  33  3  3  "  ?????  11111 

666  _  _  4644464  UU'JUt  333333  ?P???P  1,(1 

6666666664444  6444444  4  4  U  U  U  4  333333  p????? 

46646466444464646444  4444466  u  J  4  4  u  U  3333333  ?P?? 

646444  44446444444  6666546  UUdUUUO  3  3  3  3  3  3*1 

6664 _  _ —  44444444  4  6  46^66  JUUUUU/IU  3} 

666  7  7  7  7  7  V"7  ~  “  66666666  44444446  "  yuuuij  iTum  uu 

666  777777777777777777  646664446  1,444646646 

66  777777 777Vf  77777 7 7 77 77 777 7 77  66666666466  5 

6V  “7r  “  ~MfiSA55«B085"aaB 

6  9QQ99QQQQQU9Q94QQQQQQQQQQQPQQQQClQQ94U,}OOoyC1OQCOoo0Q0QQQOQQQQ0QQQQQQ 


Fig.  41a  Quantized  Contours  of  a  Distribution  Function 
Calculated  Using  the  Coefficients  of  Fig.  40 
with  E  =  8  MW/cm  and  \  =  1.06  fim.  An 
initially  constant  distribution  function  collapses 
in  energy  space  for  0.  14  psecs. 
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Fig.  42  Ionization  Rate  From  Experiment  and  Calculated  by 
Including  the  Deformation  Potential.  The  solid 
curves  include  the  full  energy  loss  for  deformation 
potentials  while  the  dotted  curve  includes  deforma¬ 
tion  potential  effects  in  the  collision  frequency  but 
not  in  the  energy  loss.  The  curves  are  again 
asymptotic  to  the  dotted  lines  drawn  so  that  r  / = 
constant.  The  computed  curves  are  for  \  =  1.06  and 
10.6  pm;  the  experiment  used  1.06  pm  radiation. 


so  that  most  electrons  caused  ionization  near  I  and  that  the  results  were 
relatively  insensitive  to  changes  in  this  cross  section. 

The  difficulty  of  increasing  the  energy  of  thermal  electrons  consti¬ 
tutes  a  low  energy  trap  which  may  affect  the  statistics  u  of  the  starting 
process  for  avalanche  breakdown,  but  these  starting  statistics  were  not 
studied. 


6*  Conclusions 


In  answer  to  the  questions  posed  in  the  Introduction,  we  conclude 
that  the  experimental  results  are  in  agreement  w’+h  calculations  based  upon 
the  hypothesis  that  the  process  causing  optical  damage  to  NaCl  can  be  des¬ 
cribed  by  a  classical  diffusion  equation.  This  agreement  may  be  considered 
very  good  considering  the  incomplete  state  of  the  theory  for  calculating  the 
collision  frequency  for  hot  electrons,  and  the  fact  that  our  treatment  ignores 

effects  due  to  the  finite  energy  quanta  in  1. 06  pm  laser  radiation  and  the  un- 

15  - 1 

certainty  in  the  electron  energy  due  to  collision  frequencies  of  10  sec 

The  experimental  results  suggest  that  the  deformation  potential  in¬ 
teraction  suggested  by  Seitz  is  the  most  important  source  of  momentum 
transfer  for  electrons  with  energy  greater  than  1  eV  in  NaCl.  It  is  not  suf¬ 
ficient  to  include  only  the  polar  interaction  described  by  Frohlich  because 

15  -1 

this  implies  collision  frequencies  below  10  sec  above  1  eV  and,  if  this 
were  the  case,  the  1.06  pm  laser  would  require  much  larger  rms  break¬ 
down  fields  than  the  10.6  pm  laser. 

As  usual,  the  Fokker- Planck  equation  is  derived  under  the  assump¬ 
tion  that  the  distribution  function  can  be  expanded  in  a  Taylor  series.  As 
discussed  more  thoroughly  in  Appendix  A;  this  assumption  is  valid  for  fields 
up  to  about  20  MV/ cm  which  is  sufficient  to  describe  the  breakdown  process 
even  for  psec  pulses.  In  addition,  we  assumed  that  the  dispersion  due  to 
energy  losses  to  the  lattice  could  be  neglected  in  computing  the  diffusion 
coefficient  D;  in  Appendix  A  we  show  that  this  approximation  is  no  longer 
ligible  when  the  field  drops  below  about  5  MV/cm.  Consequently,  the 
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ionization  rates  computed  in  Fig.  42  are  in  good  agreement  with  the  exper 
rnents  involving  15  and  300  psec  pulse  durations.  17  but  the  calculated  rate 
fall  below  the  measurements  for  4.7  and  10.  3  ns  pulses. 

The  computations  also  indicated  that  electrons  often  diffuse  from 
the  low  energies  to  the  high  energies  when  the  "average"  electron  is 
actually  losing  energy  over  a  large  part  of  the  energy  range. 
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APPENDIX  A 


RANGE  OF  VALIDITY  OF  T  HE  FOKKER-PLANCK  EQUATION 

The  Shockley  theory  of  ionization78'80*83  suggests  ionization  rates 
should  be  proportional  to  exp  -  (1/  eE I )  where  i  is  the  mean  free  path. 
Such  an  expression  can  be  conveniently  plotted  on  semilog  paper  since 
log  p  is  a  linear  function  of  1/  E.  In  fact,  the  experimental  results  fit  this 
form  reasonably  well,  although  the  Shockley  slope  is  sensitive  to  i  .  Our 
calculated  values  for  p,  which  agree  well  with  the  high  field  measurements 
in  Fig.  42,  are  more  conveniently  plotted  on  log-log  plots  to  snow  that  p 
is  asymptotically  proportional  to  E2. 

The  Shockley  argument  emphasizes  the  exceptional  electrons  which 
gain  energies  equal  to  I  without  undergoing  collision.  Actually  these 
exceptional  electrons  also  contribute  to  the  calculation72  of  the  Fokker-Planck 
coefficients,  although  they  are  inaccurately  treated  by  this  approach. 

The  contribution  to  the  ionization  rate  due  to  the  exceptional  electrons 
alone  will  be  of  the  order  of  v  exp  -  (1/  eEf )  since  each  electron  undergoes 

v  collisions  per  unit  time.  This  contribution  will  be  negligible  compared 
with  the  rates  in  Fig.  42  if 

I/SE,>>1  •  (28) 

Although  t  is  a  function  of  energy,  we  can  approximate  it  by  its 

value  at  u  =  1/2  using  the  relation  (  =  (2u/ m)1/ 2/ Then  the  criterion 
of  Eq.  (  28)  is 

I2/ e2E2£2  »  1 


or 


(29) 


2„2 

e  E 


m  v 


«  I 


This  last  expression  is  exactly  the  condition  we  have  already  given  for 
applying  the  Fokker-Planck  equation  (when  u  =  0).  Thus  the  Shockley 
effect  becomes  important  for  fields  greater  than  the  breakdown  fields  meas¬ 
ured  in  NaCl.  Even  in  other  materials,  breakdown  will  often 
occur  before  Eq.  (29)  is  violated  because  ionization  will  be  caused  at 
lower  fields  by  electrons  which  suffer  many  collisions  but,  by  chance, 
they  rebound  with  their  velocity  oriented  so  that  they  continue  to  gain 
energy  from  the  field.  The  diffusion  equation  (whether  obtained  from  the 
Legendre  expansion  or  the  Fokker-Planck  approach)  can  be  expected  to 
describe  these  electrons  accurately,  when  Eq.  (  29)  is  obeyed. 


The  diffusion  coefficient  used  in  this  paper  has  ignored  the  dis- 

PerS72?76  (du/  dt)L*  an  approximation  which  is  valid  for  a  cold  background 
gas,  *  i.e.  ,  kT  «  u.  For  finite  kT  and  weak  fields,  the  diffusion 
coefficient  D  =  -^  uA  must  become 


D  =  |  uA+P 


(30) 


where  P  is  independent  of  the  field.  When  P  is  not  zero,  G  is  given  by 


G  = 


,du.  9P 
dt  L  <Ju 


(31) 


If  an  equilibrium  function,  fQ,  is  known  for  zero  field,  and  the 
ionization  rate  is  negligible  so  that  the  current  is  zero,  then  P  must 
satisfy 


9f 

J=  Gfo'D^  =  0 


(32) 


where,  taking  fQ  to  be  the  Maxwellian,  V  u  exp  (-u/kT),  Eq.  (32)  becomes 
a  differential  equation  for  P.  .  (^-  jL)  P  -  <*■>  which  has  the 
solution 


P  =  -  u"1/  2  exp  (u/kT)  Q 


(33) 


A -2 


where 


u 

Q  =  f  u'  !  /  2  exp  -  (u'  /  kT)  (^"^l  c*u' 
o 


The  integral  Q  is  proportional  to  the  average  of  the  energy  loss 
over  a  Maxwellian  distribution  and  thus  must  equal  zero  if  the  upper  limit 
is  taken  to  be  infinity. 

As  an  example  of  the  use  of  Eq.  (33),  it  is  instructive  to  consider 
an  electron  of  energy  u  colliding  elastically  with  a  collision  frequency  v 
with  gas  molecules  of  mass  M  in  a  Maxwellian  distribution  at  tempera¬ 
ture  T.  Then  the  rate  of  energy  loss  from  the  electron  to  the  molecules 
can  be  shown  to  be 

(ht)l  =  ev(u "  3kT/2)  -  eukT§ir  ’  (34) 

3 

where  the  physical  reason  for  the  last  term  is  that,  if  is  positive, 
molecules  moving  antiparallel  to  the  electron  are  more  likely  to  cause  a 
collision  and  these  molecules  have  a  higher  energy  in  a  center  of  mass 
coordinate  system  than  molecules  moving  parallel  to  the  electron.  Here 
e  =  2mM/  (m+M)2.  Substituting  Eq.  (34)  into  Eq.  (33),  Q  is  seen  to  be 
a  perfect  integral  such  that 

P  =  evkTu  ,  (35) 


as  can  be  verified  independently.  If  we  take  J  =  0  and  a  nonzero  field, 

84 

substituting  Eqs.(30)  and  (31)  into  (16)  yields 


f  =  f (0 )  exp 


A  +  3ev  (-?! —  u) 
(2A  +  3  evkt) 


(36) 


If  much  more  energy  is  gained  from  the  field  in  a  collision  time 
then  is  lost  by  a  thermal  electron  so  that 


2  2 
e  E^ 

T 


m(v  +  w  ) 


ekT 


(37) 


A-3 


i.e.  ,  so  that  A  »  P/u,  Eq.  (36)  becomes 


u 


f  =  f (0)  u  exp 


1-15 


£v 


du 


(38) 


- — -  n 

which  has  the  form  V  u  exp  (-  g(u)/  E  )  corresponding  to  Wolff's  diffusion 
77  78 

form.  ’  If  Eq.  (37)  is  satisfied,  it  is  permissible  to  assume  that  the 
Maxwellian  describing  the  background  gas  has  a  zero  temperature,  i.  e.  , 

O  it 

it  becomes  a  delta  function  as  was  assumed  by  Wannier,  in  his  version 
of  the  diffusion  equation  which  was  intended  for  sufficiently  high  fields. 


When  ionization  is  occurring  it  is  no  longer  strictly  correct  to  set 
J=0  since  the  ionization  rate  is  proportional  to  J(I),  and  Eq.  (36)  must 
be  considered  an  upper  bound  as  in  Eq.  (17).  However,  it  is  probably 
true  that  Eq.  (38)  is  an  approximation  to  the  distribution  function  so  that 
the  ionization  rate  is  proportional  to  Gf  -  D  i.  e.  ,  proportional  to  E2 
as  the  numerical  solutions  show.  Of  course,  the  exact  values  for  the 
ionization  rate  cannot  be  obtained  from  a  solution  which  assumes  that 
J  =  0,  which  prevents  us  from  using  the  elegant  formalism  of  Keldysh80 
to  obtain  accurate  results. 


If  we  take  T  =  300  K,  €  =  10  ^  and  2  +  u2  =  1080,  the  inequality  in 
Eq.  (37)  becomes  an  equality  at  .05  MV/cm. 


It  is  tempting  to  find  Q  for  NaCl  by  substituting  (du/  dt)j^  into 
Eq.  (33).  However,  if  we  substitute  (du/  dt)L  =  d1  u-1'  2fn  d2  u  into 
Eq.  (33),  we  obtain 

u/kT 

Q  =  d^kT  in  d2kT[  1  -  exp  (-  u/ kt)]  +  \  e“x  in  x  dx 

o 


(39) 


which  does  not  go  to  zero  for  largeu  ifweuse  thevalue  of  d2  in  Eq.  (18).  The  reason 
is  that  the  derivation  of  Eq.  (18)  is  only  accurate  for  energies  greater 
than  thermal  energies,  and  (du/  dt)j^,  is  negative  below  thermal  energy 


A-4 


in  order  to  balance  energy  gains  and  losses  in  the  equilibrium  distri¬ 
bution.  We  make  an  arbitrary  fixup  by  choosing  d  so  that 

oc 

£nd9kT  =  -  ^  e  x  fn  x  dx  , 
o 


when  we  calculate  P.  Then 


l 


P  = 


djkT  u"1/ 2 


-x 

e 


fnxdx 


oc 

f  -x 

\  e  in  x  dx 
o 


4 

where  y  -  u/kT,  which  can  be  numerically  integrated  to  obtain  the  values 
in  Table  II.  The  inequality  A  »  3P/  2u  requires  E  to  be  greater  than 
0.  15  MV/  cm  assuming  vj  (v2  +  w2)  =  lO-15,  u  =  1/ 2  and  P  =  0.  1  eV2/sec. 


We  have  added  P  to  the  diffusion  term  in  the  Fokker-Planck  equa¬ 
tion,  with  the  same  value  for  (du/  dt)L  that  we  used  previously.  The  result 
for  a.  =  10.  6  pm  and  E  =  7  MV/  cm  in  Fig.  42  which  had  been  8.8  XlO^sec'1 

was  increased  by  8.2  percent.  The  increasein  q  is  corresnondingly  greater 
for  fields  below  7  MV/ cm. 


TABLE  II 

u(eV)  .1  .5  l.i  2.1 

7^2  ; - 

^psec*  1  -192  -144  -H8  :  .097 


3.1  4.1  6.1  I  8.1 

.086  .079  .069  .063 


The  values  of  P  given  here  are  only  a  rough  approximation  since 
they  depend  on  a  model  of  the  energy  loss  to  the  lattice  which  we  know  to 
be  inexact  near  thermal  energies. 


Our  arguments  are  not  intended  to  deny  that  the  low  field  data  for 
n  may  be  approximately  described  by  a  curve  proportional  to  exp  (-constant/  E) 
since  our  experimental  results  lie  close  to  such  a  curve,  but  only  that  the 
low  field  behaviour  may  be  explored  by  means  of  the  Fokker-Planck  equa¬ 
tion  without  undertaking  a  special  treatment  for  those  exceptional  electrons 
which  do  not  suffer  any  collisions  before  they  reach  the  ionization  energy. 


A -5 


APPEND  LX  B 


PUBLICATIONS  AND  PRESENTATIONS  ON  VISIBLE  AND  INFRARED 
LASER-INDUCED  DAMAGE  TO  TRANSPARENT  MATERIALS 


1.  I).  W.  Fradin  and  M.  Bass,  "Electron  Avalanche  Breakdown  at  Optical 
Frequencies,  "  Winter  Meeting  of  the  American  Physical  Society, 
paper  GE3,  New  York,  February  1973. 

2.  D.  W.  Fradin  and  M.  Bass,  "A  Comparison  of  Laser-Induced  Surface 
and  BuIk  Damage,  "  Appl.  Phvs.  Lett.  22_,  157  (1973). 

3.  D.  W  .  Fradin  and  M.  Rass,  1  Electron  Avalanche  Breakdown  Induced  by 
Ruby  Laser  Light,  "  Appl.  Phys.  Lett.  22_,  206  (1973). 

4.  D.  W.  Fradin,  E.  Yablonovitch,  and  M.  Bass,  "Confirmation  of  an 
Electron  Avalanche  Causing  Laser-Induced  Bulk  Damage  at  1.06  Microns," 
Appl.  Optics  V2,  700  (1973). 

5.  D.  W.  Fradin,  "The  Relationship  Between  Self-Focusing  and  Optical  Bulk 
Damage  and  the  Measurement  of  Self-Focusing  Parameters,"  Twelfth 
Symposium  on  Electron,  Ion,  and  Laser  Beam  Technology ,  MIT, 
Cambridge,  Mass.  (1973). 

6.  D.  W.  Fradin  and  M.  Bass,  "Studies  of  Intrinsic  Optical  Breakdwon," 
Proc.  5th  ASTM-NBS  Symposium  on  Laser  Damage,  Boulder,  Colorado, 
May  15-16,  1973  (NBS  Spec .  Publ.  387). 

7.  M.  Bass  and  D.  W.  Fradin,  "Statistics  in  Surface  and  Bulk  Breakdown 
and  Identification  of  Intrinsic  Mechanisms,"  IEEE  J.  Quant.  Elect. 

QE-4,  890  (1973). 

8.  M.  Bass,  "Laser  Damage  to  Transparent  Media,  "  Invited  Paper  19.6, 

1973  Conference  on  Laser  Engineering  and  Applications,  Digest  of 
Technical  Papers,  Washington,  D.  C.  ,  May  30  -  June  1,  1973. 

9.  D.  W.  Fradin,  'The  Measurement  of  Self-Focusing  Parameters  Using 
Intrinsic  Optical  Damage,  "  IEEE  J.  Quant.  Elect.  QE-9,  954  (1973). 

10.  D.  W.  Fradin,  M.  Bass,  D.  P.  Bua,  and  C.  A.  Christian,  "Optical 
Damage  in  Laser  Materials,  "  Proc.  of  Air  Force  Conference  on  High- 
Power  IR  Laser  Window  Materials,  Hyannis,  Mass.  ,  November  12-14, 
1973. 

11.  D.  W.  Fradin,  N.  Bloembergen,  and  J.  P.  Letellier,  "Dependence  of 
Laser-Induced  Breakdown  Field  Strength  on  PUlse  Duration,"  Appl. 

Phys.  Lett.  22_,  635  (1973). 

12.  D.  W.  Fradin  and  M.  Bass,  "Effects  of  Lattice  Disorder  on  the  Intrinsic 

Optical  Damage  Fields  of  Solids ,  "  Appl.  Phys.  Lett.  23,  604(1973). 


B  - 1 


1.-1.  H avid  Milam,  R.  A.  Bradbury,  and  M.  Bass,  "Laser  Damage  Threshold 
lor  Dielectric  C  oatings  a0  Determined  by  Inclusions,"  Appl.  Phvs.  I  ett 
££,  6d4  (1973). 

14.  D.  W  .  1-rad  in  and  D.  P.  Bua,  "Laser-Induced  Damage  in  ZnSe ,  "  submitted 
for  publication  (Feb.  1974). 


